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SUMMARY 


A program called FLEXWAL is presented. The program was 
written to predict the wall modifications necessary to remove the 
upper and lower wall interference effects in solid, flexible, 
adaptive-wall wind tunnels. The procedure used is the iterative 
method of NASA TP-2081 and is valid for both subsonic and transonic 
test conditions. Instructions on program usage, functional state- 
ments of program routines, a program listing, and a sample case are 
given. 
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INTRODUCTION 


The computer program FLEXWAL described in this report is used 
to calculate the required upper and lower boundary modifications 
necessary for removing wind-tunnel wall interference effects in 
solid flexible-wall wind tunnels. The method used is the iterative 
technique of reference 1 and requires experimentally obtained 
information as input. The program is completely self-contained and 
is written in a modular form so that changes to the program for 
specific applications or wind-tunnel facilities can be made with 
ease. 

Both analytical and experimental applications of the program 
at subsonic and transonic speeds have been previously presented 
along with the convergence properties of the method (See ref. 1.). 
This report gives a general description of the method, instructions 
for using the program, sample input, functional statements of 
routines used, and a listing of the computer code. 
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GENERAL DESCRIPTION 


When testing in a wind tunnel, the flow field about the model 
is constrained to something different than that which would be 
present in free air. The principal difference is referred to as 
the wall interference and is generally taken as a perturbation 
about the free-air case. Classical methods for determining the 
wind-tunnel wall interference make mathematical approximations of 
the perturbation from which correction formulae are derived. The 
present method is that of reference 1 and uses the Cauchy integral 
formula to analytically extend the real flow field in the wind tun- 
nel to infinity by solving for an imaginary flow exterior to the 
wind tunnel. These two flows are then iteratively coupled at the 
boundary with continuity being established when the measured data 
match the results of calculated data at the wall. If there is a 
mismatch then the method predicts a correction to the wind-tunnel 
wall, new data are measured, and the continuity is again checked. 

Functional statements of the routines used in the FLEXWAL 
program are given in Appendix A and a listing of the program is 
contained in Appendix B. The program is completely self-contained 
and modular in construction. All input data are read into subrou- 
tine INPT and an echo listing of the input is given in the output 
listing. Subroutine SCALE uses the Goethert scaling laws to remove 
the effects of compressibility at the test Mach number prior to 
application of the correction technique and then, after completion 
of the correction calculations, SCALE is called again so as to 
include the compressibility at the desired Mach number MSET. MSET 
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has been included for consistency during the iterative process 
because often the desired Mach number is slightly different than 
that recorded during the test due to experimental error. Subrou- 
tine SMOOTH is a simple least squares smoothing subroutine which 
may or may not be needed depending on the quality of the experi- 
mentally obtained data. For each set of five adjacent pressures, a 
least squares parabola is determined and the center pressure is 
calculated from the least squares solution. When all pressures 
have been smoothed, a single pass has been completed and the pro- 
cess will be repeated as specified by the IS parameter in the 
input. The subroutine WALLDTA determines values of the perturba- 
tion velocities and sets up an integration contour as specified in 
reference 1 . A spline through the wall jack coordinates allows the 
wallposition to be determined at intermediate stations and by 
taking the derivative of the spline curve, the wall slopes (i.e. 
vertical velocities) may also be obtained along the walls. In the 
present form of the computer code, the vertical velocity components 
on the upstream and downstream legs of the contour have been set 
equal to the average of those on the upper and lower walls at the 
same station. This can easily be changed if the appropriate data 
are measured or better approximated. 

The actual application of the method is in subroutine VJACK. 

It is here that the contour integration is performed for the 
velocities at the wall jack stations. It is important to note that 
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the most upstream and downstream wall jacks do not actually have to 
lie on the data contour. The program will determine whether the 
velocity is being calculated at a field point or a point on the 
contour and will adjust the solution accordingly. The difference 
is a factor of 2 which can be seen by referring to reference 1 . 
Subroutine DELYJ calculates the predicted wall change by integrat- 
ing slopes determined from the new jack velocities along each 
wall. Convergence of the method is declared when these changes 
become sufficiently small or when some other criteria such as 
negligible variations in the airfoil force coefficients is met. 
Subroutine NEWYJ prints the results of the program. All other 
subroutines act in support of the previously noted routines. 

PROGRAM USAGE 

The program is written in the FORTRAN IV computer language and 
requires 42000 words of execution field on a CDC CYBER 173 com- 
puter. Typical program run times are about 2 seconds. Appendix C 
gives a sample case. The output listing gives an echo listing of 
the input data required to generate the results. A user's guide 
containing the input format, a glossary of all terras used in the 
input and output, and program restrictions are presented in the 
computer program and are found at the top of the program listing in 
Appendix B. 
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APPENDIX A 


SUMMARY OF PROGRAM ROUTINES AND THEIR FUNCTIONS 



ROUTINE 

NAME 

FLEXWAL 

CURV'D 

CURVI 

CURV2 

IB DRY 
INTRVL 

ISTRIF 

RMACHNO 

RMS 

SWITCH 


EXTERNAL 

REFERENCES 


AVGDY 

DELYJ 

INITIAL 

INPT 

NEWYJ 

RMACHNO 

SCALE 

SMOOTH 

UINFIN 

VJACK 

WALLDTA 

INTRVL 

SNHCHS 


INTRVL 

SNHCHS 


INTRVL 

SNHCHS 


FUNCTIONAL STATEMENT 


This is the main or executive routine of 
a FORTRAN IV computer program for calcu- 
lating the jack settings of the upper and 
lower walls of a two-dimensional solid 
flexible wall wind tunnel. 


This function subprogram differentiates 
a curve at a given point using a spline 
under tension. The subroutine CURVI 
should be called earlier to determine 
certain necessary parameters. 

This function subprogram integrates a 
curve specified by a spline under tension 
between two given limits. The subroutine 
HMCURV1 should be called earlier to 
determine certain necessary parameters. 

This function subprogram interpolates a 
curve at a given point using a spline 
under tension. The subroutine HMCURV1 
should be called earlier to determine 
certain necessary parameters. 

This function subprogram determines 
whether a point is on the data surface 
or a field point. 

This function subprogram determines the 
index of the interval (determined by a 
given increasing sequence) in which a 
given value lies. 

This function subprogram determines 
whether the strip being considered is 
a normal or singular strip. 

This function subprogram calculates the 
test Mach numbers. 

This subroutine function determines the 
RMS value of array A. 

This subroutine reorders the indices of 
array x from high to low. 
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ROUTINE 

NAME 

UCAL 

UINFIN 

AVGDYJ 

DELYJ 

HMCURV1 

INITIAL 

INPT 

LSQSMTH 

NEWYJ 

SCALE 

SING 

SMOOTH 

SNHCSH 

STRIP 


EXTERNAL 

REFERENCES 


CURVI 

HMCURV1 

RMS 


RMS 


LSQSMTH 


FUNCTIONAL STATEMENT 


This function subprogram calculates the 
U component of velocity. 

This function subprogram calculates the 
free-stream velocity. 

This subroutine averages the top and 
bottom wall values of DYJ. 

This subroutine integrates the delta 
theta values to obtain new wall 
location predictions. 

This subroutine fits a cubic or tension 
spline through a set of Y vs. X values. 

This subroutine presets certain arrays to 
zero. 

This subroutine handles the input of the 
raw data as taken from the wind-tunnel 
test . 

This subroutine provides a least squares 
smoothing for either a straight line or 
a parabolic fit of the input data. 

This subroutine writes the predicted 
results onto TAPE 6. 

This subroutine scales the input using 
the Goethert scaling rules. 

This subroutine evaluates the velocity 
component along the singular strip. 

This subroutine controls the smoothing 
of the wind-tunnel data. 

This subroutine returns approximations 
to: SINHM(X) = SINH(X) - X 

COSHM(X) - COSH(X) - 1 

and 

COSHMM(X) = COSH(X) - 1 - X*X/2 
with relative error less than 3.24E-14. 

This subroutine calculates the jack 
velocity component resulting from some 
strip excluding that containing the 
singularity. 
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ROUTINE 

NAME 

VJACK 

WALLDTA 


EXTERNAL 

REFERENCES 


FUNCTIONAL STATEMENT 


IBDRY 

ISTRIP 

STRIP 


This subroutine integrates the velocity 
integral to obtain the jack velocities. 


CURVD 

CURV2 

HMCURV1 

SWITCH 

UCAL 


This subroutine takes input data and 
calculates velocities and wall locations 
then fills appropriate array locations. 


10 



APPENDIX B 


PROGRAM LISTING 



c*«» ( »«»««»«««„„„„ 4 „„ 9t5 „ # „„„„„ s „ JM4}{MS4MM4MMtM5tM 

c * 


c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 


FLEXWAL IS a FORTRAN IV COMPUTER PROGRAM FOR CALCULATING THE JACK 
SFTTIN6S OF THE UPPER AND LOWER WALLS OF A TWO-DIMENSIONAL SOLID 
FLEXIBLE-WALL WIND TUNNEL AS PREDICTED BY THE METHOD OF NASA TP-2081. 
THF PROGRAM WAS WRITTEN BY JOEL L. EVERHART OF NASA LANGLEY RESEARCH 
C t MTER • 


REQUIRED INPUT IS READ INTO THE PROGRAM FROM TAPES AND OUTPUT IS 
RETURNED ON TAPES. 


# 

•a- 

■fr 

■& 

* 

a- 


C *********** ************ ***«***#«.** 
C* 


■frtftf’&tttt-H-* *> ************ ■t:-*^************^^.^. 


c* 

c* 

INPUT: 

c* 

c* 

READ ORDER 

c* 

1 

c* 

2 

c* 

3 

c* 

c* 

4 

c* 

5 

c* 

6 

c* 

7 

c* 

8 

c* 

c* 

9 


DESCRIPTION 

HFAD(I) ,1=1,8 
TEST * RUN* PO I NT 

PNf* ALPHA, TEMP, PT, PINF, MSET 

N'JT , NJB, NTW, NDS , NBW , NUS, NRL A X , 

<IS(I) ,1=1,4) , I A VG , I PRT 

I OPT ( I ) ,1=1,4 

EC ( I ) ,1=1,4 

RLAX(I) ,1=1, NRL AX 

(XJ(I) , Y J ( I ) ,1=1, NJ) 

(XD(I) ,YD(I) ,UD(I) ,1=1, NO) 


FORMAT 

BA 1 0 
315 

6F10.0 

1315 

415 

4F10.0 
8F10.0 
2F10.0 
3F 1 0 • 0 


« * * ■» * ** * * 
«■ 
* 

* 

* 

* 

*> 

*• 

* 

* 

* 

* 


C* 


* 

***»»## 


C* DEFINITION OF TFPMSJ 
C* 


C* 

ALPHA 

c* 

DUJ 

C* 

DVJ 

C* 

DYJ 

c* 

EC 

C* 


C* 


c* 


c* 


c* 


c* 


c* 

HEAD 

c* 

I AVG 

c* 



MODEL ANGLF OF ATTACK 

PREDICTED CHANGE IN AXIAL JACK PERTURBATION VELOCITY 
PREDICTED CHANGE IN VERTICAL JACK PERTURBATION VELOCITY 
PREDICTED CHANGE IN WALL POSITION 
END CONDITIONS ON THE SPLINE THROUGH THE JACK 
LOCATIONS. THE SPLINE IS USED TO DETERMINE WALL POSITION 
BETWEEN JACKS AMD WALL SLOPES. 

EC ( 1 ) » UPSTREAM.UPPER wall 
EC ( 2) * DOWNSTREAM, UPPER WALL 
EC 1 3 ) , UPSTREAM, LOWER WALL 
EC ( 4 ) , DOWNSTREAM, LOWER WALL 
DESCRIPTOR CARD OF UP TO 80 CHARACTERS. 

WALL JACK AVERAGING PARAMETER. 

AT A GTVFM XJ STATION THF UPPFP AND I OWFR YJ POSTTTON ARF 


* 

* 

# 

«& 

* 

* 

* 

* 

* 

* 



1 

1 

c* 



AyERAGFD giving SYMMETRICAL results 

* 


c* 



iavg=o» averaging on 



c* 



1* AVERAGING off 

* 


c* 

I OPT 

* 

SPLINE FNO CONDITION OPTION PARAMETER 

# 


c* 



ICPT(l), UPSTREAM, UPPER WALL 

a 


c* 



I OPT ( 2 ) , DOWNSTREAM, UPPER WALL 

a 


C# 



I OPT ( 3 ) , UPSTREAM, LOWER WALL 

# 


C* 



I OPT ( 4 ) , DOWNSTREAM, LOWER WALL 

a* 


C* 



I OPT { I ) =0 , EC UNKNOWN 

a 


c* 



=1, EC=D(Y)/DX 

a 


c* 



=?, EC=D2(Y)/0X2 

a 


c* 

IPRT 

- 

PRINT CONTROL PARAMETER. 

a 


c* 



IPRT=0» MAXIMUM OUTPUT 

a 


c* 



=1. NORMAL OUTPUT 

a 


c* 

IS 

_ 

LEAST SQUARES SMOOTHING PARAMETER FOR INPUT DATA ON 

a 


c* 



FACH WAIL OF CONTOUR. I S ( I ) IS THE NUMBER OF SMOOTHING 

a 


c* 



PASSFS MADE ON EACH WALL STARTING CLOCKWISE ON UPPER WALL. 

a 


c# 



IS(I)=0 IMPLIFS NO SMOOTHING 

a 


c* 

MSET 

• 

DESIRED MACH NUMBER TO WHICH RESULTS ARE SCALED. 

a 


c* 

HD 

- 

TOTAL NUMBER OF PRESSURE MEASUREMENTS (CALCULATED 

a 


c* 



INTERNALLY TO PROGRAM) 

a 


c* 

NRVI 

— 

NUMBER OF PRESSURE MEASUREMENTS ON BOTTOM WALL (LOWER LEG 

a 


c* 



OF CONTOUR) 

a 


c* 

NiDS 

_ 

NUMBER OF PRESSURE MEASUREMENTS ON DOWNSTREAM LEG OF CONTOUR 

a 


c* 

NT W 

- 

NUMBER OF PRESSURE MEASUREMENTS ON TOP WALL (UPPER LEG OF 

a 


c* 



CONTOUR) 

a 


c* 

NIJS 

_ 

NIJMRRR OF PRESSURE MEASUREMENTS ON UPSTREAM LEG OF CONTOUR 

a 


c* 

NJ 

— 

TOTAL NUMBER OF JACKS (CALCULATED INTERNALLY TO PROGRAM) 

a 


c* 

N JR 

- 

NUMBER OF JACKS ON BOTTOM WALL 

a 


c* 

NJT 

- 

NUMBER OF JACKS ON TOP WALL 

a 


c* 

NRLAX 

• 

NUMBER OF WALL RELAXATION FACTORS 

a 


c* 

PINF 

- 

FREESTREAM STATIC PRESSURE, RSI OR PSF 

a 


c* 

POINT 

- 

DATA IDENTIFIER 

a 


c* 

PT 


FRFESTPFAM TOTAL PRESSURE, PSI OR f>SF 

a 


c* 

RLAX 

— 

VALUE OF WALL CHANGE RELAXATION FACTORS FOP FINAL RESULTS 

a 


c* 

PMS 

— 

ROOT-MEAN-SQUARE ERROR 

a 


c* 

PM 

~ 

TEST RFYNOLDS NUMBER 

a 


c* 

PUN 

- 

data identifier 

a 


c* 

TEMP 

- 

FREESTPEAM TOTAL TEMPERATURE, DEG. RANKINE 

a 


c* 

TEST 

« 

DATA IDENTIFIER 

a 


c» 

THETA 

- 

flow angle at wall 

a 


c* 

THETAM 

— 

PREDICTED FLOW ANGLE AT WALL 

a 


c* 

UD 

• 

ON INPUT, MEASURED WALL PRESSURE COEFFICIENTS. ON OUTPUT, 

a 


c* 



CALCULATED AXIAL VELOCITY. 

a 


c* 

UJ 

- 

AXIAL VFLOCITY AT JACK 

a 


c* 

IJJN 

- 

PREDICTED AXIAL JACK VELOCITY 

a 

LJ 

c* 

VJ 

— 

VERTICAL VELOCITY AT JACK 

a 

W 

c* 

VJN 

• 

PREDICTED VERTICAL VELOCITY AT JACK 

a 


c* 

XD 

- 

AXIAL ST AT I ON OF WALL ORIFICES 

a 


c* 

XJ 

— 

AXIAL STATION 0^ WALL JACKS 

a 


c* 

YO 

- 

MOMTNAI VFRTTCAI 1 or. AT TON OF IIPPKR AND 1 owfr WAI 1 . ORTFTCFS. 

a 



ACTUAL VERTICAL LOCATION 1 OF UPSTREAM AND DOWNSTREAM ORIFICES. 
ACTUAL UPPER and LOWER wall ORDINATES ARE CALCULATED FROM A 
SPLINE THROUGH THF JACK ORDINATES. 

VERTICAL LOCATION OF WALL JACKS. 

PPF 0 1 CTED VERTICAL LOCATION OF WALL JACKS. 


C* 

c* 
c* 

C* YJ 

C* YJN 

C* 


» 

* 

» 

* 

« 


RESTRICTIONS: 


?. 

3. 

4 . 


WALL DATA ARE INPUT CLOCKWISE STARTING AT THE MOST UPSTREAM UPPER 
WALL ORIFICE. 

THF JCCK ORDINATES AOF INPUT FROM UPSTREAM TO DOWNSTREAM URPFP WALL 
AND FROM UPSTREAM TO DOWNSTREAM LOWER WALL. 

USE CONSISTENT PRESSURE UNITS. 

SMOOTHING OF UPSTRFAM AND DOWNSTREAM DATA IS NOT INCLUDED. 

HOWFVFR, NOTE IS MADE IN THF LISTING OF SUBROUTINE. SMOOTH WHERE 
IT SHOULD BE INCLUDED IF DESIRED. 

THE FLOW ANGLES (IE. VD/UD) ON THE UPSTREAM AND PtOWNSTRE am 
PORTIONS OF THE DATA CONTOURS ARE PRESENTLY SET F DUAL TO THE 
AVERAGE OF THF SLOPES OF THE UPPER AND LOWER WALLS AT THAT POSITION. 


* 

•» 

* 

* 

* 


c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 

C *************************************************** ***************************** 
->#*##£• ***#*##*******************************^^'^^*'^^^****^^^*’'^*^^*^^ 

c 

c 

c 

c 

c 

c 

PROGRAM FLEXteAL ( INPUT * OUTPUT t TAPES* TAPES) 

REAL M t M 5? E T 

INTF6ER TEST * RUM f POINT 

COMMON/ J A CK/NJ.NJT.NJP ♦ X J ( 1 00 ) » YJ ( 1 00 ) »UJ ( 1 00) .VJ(100) . 

SU JN (100) .VJN(IOO) .YJN(IOO) .DYJ(IOO) 

COMMON/P ARAM/HE AD (8) » TEST * RUN .POINT tM ♦ RN . ALPHA * TEMP. PL AX ( 10) ♦ 

"SKIRL AX .PT.PINF .UIMF .IAVG»IS(4) »IOPT(4) .EC(4) . GAMMA . MSF T 
COMMON /CO NT OUP /ND . NTW . NDS ♦ NBW . NUS .XD(100) * YD ( 1 00 ) .UD < 100) . VU < 100 ) 

DIMENSION WK1 ( 1000) »WK2 ( 100 ) 

«AMMA=1 .4 
OASP= 1 7 1 A . 0 


NWK1=1000 

nwk?=i oo 

CALL INITIAL (NWK2.WK2) 

CALL INPT(IPPT) 

MsRMACHMO (GAMMA.pt. pi NF ) 
UINF=!JINFTN (GAMMA .GASP. TEMP «M) 
WRITE (G.Q6) M.UINF 
IF ( IPRT.FQ.O) WRITE (6.RR) 


25 

96 

97 

98 

99 


CALL SCALE (NJ*XJ* YJ*DYJ*UJ* VJ*M*MSET * 0 * 1 * IPBT) 
IF(IPRT.EO.O) WRITE (6*98) 

CALL SCALE (ND*XP*Y0*WK2*UD* VD»M,MSET» 0*2* IPRT ) 
CALL SMOOTH ( 50 0 * VK 1 (1> » WK 1 (501) »IS*IPRT) 

CALL WALLDTA(NWK1*WK1*IPRT> 

CALL VJACK(IOPT) 

CALL DELYJ (NWK2*WK2»NWK1 *WK1 *UINF ) 

IF(IPRT.FO.O) WRITE (6*97) 

CALL SCALE (NJ*XJ* YJ»DYJ*UJN» VJN*M*MSET* 1 *3* IPRT) 
IF(IAVG.EQ.O)CALL AVGDYJ(MJ*XJ»DYJ) 

00 25 I=1,NPLAX 

CALL NEWYJ(I.NWK2,WK2) 

CONTINUE 

STOP 

FORMAT (///* M*INF=*G12.6* U I NF=*G 1 2 . 6// ) 

FORMAT («1 m //* CALCULATED JACK DATA*) 

FORMAT ("1"//* INPUT WALL DATA*) 

FORMAT ("1"//* INPUT JACK DATA*) 

END 


cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


c 

c 

c 


ccc 
c 
c 
c 


THIS FUNCTION CALCULATES THE TEST MACH NUMBER 
FUNCTION RMACHNO ( GAMMA *PT »PINF) 

REWl w 2 

PR=(PT/PINE)**( (GAMMA-1 . 0 ) /GAMMA ) 

(PR-1 . 0) *2.0/ (GAMMA- 1.0) 

IF (M2.LT.0.0) M2=0.0 
pmACHNO=SORT (M2) 

RETURN 
FNO 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

THIS FUNCTION CALCULATES THE FREESTREAM VELOCITY 


FUNCTION UINFIN( GAMMA *GASR« TEMP *M) 

REAL M 

TINF=TEMP/ (1.0* (GAMMA-1 . 0 ) *M*M/2. 0 ) 

A2=GAMMA*GASP*TINF 
A=SORT (A2) 

UINFI M=A*M 

RETURN 

END 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c THIS SUBROUTINE PRESETS CERTAIN ARRAYS TO ZERO 


SUPROIJTTNF INITIAL (NWK2«WK2) 



o o o o n 


COMMON/JACK/N J»NJT *NJR » X J ( 100 ) » Y J ( 1 00 ) ♦UJ(100) , V J ( 1 0 0 ) » 

SUJNdOO) » V JN ( 1 0 0 ) » Y JN ( 1 0 0 ) ,DYJ(100) 

COMMON/CON TOUR /NO »NTW,NOS,NRW, NI.IS » X D ( 1 0 0 ) , YD(100)*UD(100),VD<100) 
DIMENSION WK2(NWK2) 


DO 10 

I 

= 1*100 

XJ(I) 

= 0 

.0 

Y J ( I ) 

= 0 

• 0 

UJ(I) 

= 0 

.0 

VJ(I) 

= 0 

.0 

U JN ( I 

) = 

0.0 

V JN ( I 

) = 

0.0 

YJN ( I 

>' = 

0.0 

DYJ(I 

) = 

0.0 

XD ( I ) 

= 0 

• 0 

YD ( I ) 

= 0 

• 0 

IJD ( I ) 

= 0 

.0 

VD ( I ) 

= 0 

.0 

WK2 ( I 

) = 

0.0 


10 CONTINUE 
RETURN 
END 
C 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


THIS SUBROUTINE HANDLES THE INPUT OF THE RAW DATA AS TAKEN FROM 
THE WIND TUNNEL TEST 

SUBROUTINE INPT(IPRT) 

REAL M,vsET 

INTFGER TEST»RUN»POINT 

COMMON/JACK/NJ,NJT»NJB,XJ<100) ,YJ(100) » U JO (10 0) * V JO (100) . 

SUJ(IOO) ,VJ(100) ,YJN(100) 

COMMON/P ARAM/HEAD (8) * TEST * RUN »POINT»M*RN* ALPHA, TEMP, RLAX (10) , 
SNRLAX »PT * PINF * UI NF * I A VG , IS (A) , I0PT(4) ,EC(4) , GAMMA, MSET 
COMMON/CONTOUP/ND*NTW,NDS»NBW,NUS»XD(100> »YD ( 100) *UD( 100) *VD ( 100) 
REWIND 5 

READ (5» 290) (HEAD ( I ) ,1 = 1 ,8) 

READ (5.291) TEST* RUN* POINT 

READ (5» 292) RN* ALPHA* TEMP, PT»PINF*MSET 

READ (5,291) N JT , N JB » NTW » NDS » NBW ♦ NUS ♦ NRL AX » (IS(I)*I=1«4) , IAVG,IPRT 

NJ=NJT ♦ NJR 

NO=NTW*NDS+NBW+NUS 

READ (5, 291) (IOPT(I) ,1=1,4) 

READ (5*29?) (EC ( I ) , 1=1 *4 ) 

READ (5, ?92) (RLAX(I) ,I=1,NRLAX) 

READ (5,294) (X J ( I ) ,YJ ( I ) , 1 = 1 ,NJ) 

RFAD (5*299) (XD ( I ) , YD ( I ) »UD ( I ) » 1=1 ,ND) 



o n o 


C**e**ECHO CHECK OF INPUT DATA 
C 

WRITE (6*2 95) 

WRITE (6*296) (HFAD(I) »I=1*«> 

WRITE (6.297) TEST * RUN * POINT 

WRITE (6*298) KN» ALPHA* TEMP* PT » PI NF ♦ MSE T 

WRITE (6* ?97) NJT »NJB »NTW*NDS*NBW*NUS*NRLAX » (I$(I)»I=1»4)»IAVG» IPRT 
WRITE (6*?97) ( I OPT ( I ) *1=1.4) 

WRITE (6*298) (EC (I). 1 = 1*4) 

WRITE (6*298) (RLAXCI) . I=1*NRLAX) 

WRITE (6*288) (XJ(I) *YJ(I> »I = 1*NJ) 

WRITE (6*287) (XD ( I ) »YO ( I ) »UD ( I ) » 1=1 »ND) 

RETURN 

287 FORMAT (3X«3F10»5) 

288 FORMAT (3X.2F10. 5) 

290 FORMAT (8*10) 

291 FOPMAT (1615) 

292 FORMAT (8F10*0) 

294 FORMAT (2F10»0) 

295 FORMAT ( lHl*//»* INPUT DATA LISTING **/) 

296 FORMAT ( 3X »8A 1 0 ) 

297 FORMAT (3X»16I5) 

298 FORMAT (3X»8G11#5) 

299 FORMAT (3E10»0) 

END 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


THIS SUBROUTINE CONTROLS SMOOTHING OF THE WIND TUNNEL DATA 
SUBROUTINE SMOOTH (NWK »WK 1 * WK2* IS* IPRT ) 

COMMON/CONTOUR /ND* NT W* NOS * NBW * NUS * XD (100) *YD(100) *UD (100) *VD< 100) 
DIMENSION IS (4) ,WK1 (NWK ) *WK2 (NWK ) »DUM (10) 

N1=IS(1) 

N2=IS (2) 

N3=IS (3) 

N4=IS(4) 

C**»**SMOOTH TOP WALL DATA 
IF (Ml ,FO» 0) GO TO 26 
DO 5 I=1*NTW 
WK1 ( I ) =XD (I ) 

WK2 ( I ) =HD ( I ) 

5 CONTINUF 
DO 15 I =1 *N 1 

M CALL LSOSMTH (5»WK1 ( 1 ) »WK2 ( 1 ) *5*DIJM» 1 »WK1 ( 1 ) *WKS* 0 ) 

"J WK2 ( 1 ) = WKS 

CALL LSOSMTH ( 5 » WK 1 ( 1 ) »WK2(1)*5»DUM*1»WK1(2)*WKS*0) 

WK2 (2)=WKS 
NS=NT W-2 
DO 10 J=3*MS 
T. !=.!-? 



call LS0SMTH(5,W«l (IJ) ♦WK2(IJ> *5*DUM,l,WKl (J) »WKS»0) 
WK2(J)=WKS 
10 CONTINUE 

I J=NTW-4 

KlS = (ViS*l 

CALL LSOSMTH <5» WK 1 (IJ) ♦WK2(IJ) ,5*DUM*l»WKl ( NJS ) »WKS»0) 

WK?(NS)=WKS 

MS=NS*1 

CALL LSOSMTH (5^'Kl (IJ) ♦UK2(IJ) »5»DUM,l,WKl (NS) *WKS»0) 
HK?(NS)=WKS 
15 CONTINUE 

IF(IPPT.FQ.l)GO TO 26 
00 20 1=1 ♦NTW 

VP I TF( 6. 1001X0 (I) » 1)0(1) *WK2(I) 

X0(I)=VK1 (I) 

II0(I)=VK?(I) 

20 CONTI MUF 

C*»***SMO0TH DOWNSTREAM OATA 
26 IF < N2»FO. 0 ) 60 TO 51 
C***»*SMO0TH BOTTOM HALL DATA 
51 IF(M3.FO.O)KO TO 76 
00 55 1=1 ♦ NSW 
I I=NTP+NDS+NPW*1-I 
WK1 (I)sXO(II) 

WK? ( I ) =UD (II) 

55 CONTINUE 

00 65 1=1 »N3 

CALL LSOSMTH (S,WK1 (1) ,WK2( 1 ) «5»DUM» 1»WK1 ( 1 ) «WKS*0) 
WK2(1)=WKS 

CALL LSOSMTH (5.WK1 ( l ) ♦WK? ( l ) «5»DUM, 1 » MKl ( 2) * WKS* 0 ) 

WK? (?) =‘*'KS 
NS=MBW-? 
no 50 J=3 ♦NS 
IJ=J-2 

CALL LSOSMTH (5*VK1 ( I J) »UK2 ( I J) ♦5*0UW» 1 ♦WK 1 ( J) ♦WKS^O) 
HK2(J)=WKS 
50 CONTINUE 
IJ=N6W-4 
NS=NS+1 

CALL LSOSMTH (5^WK1 ( I J) ♦WK? < I J) ♦ S.DUM* NWKi (US) ♦WKS^O) 
WK2 (NS) =WKS 
MS=NS* 1 

CALL LSOSMTH (5.WK1 (IJ) »WK2(IJ> ♦ 5 ♦ DUM* 1 ♦ WK 1 (NS) ♦WKSfO) 
WK2(NS)=WKS 
65 CONTINUE 

IF (IPRT.EO.l)GO TO 76 
00 70 1=1 ♦ NOW 
II=NTW+N0S+NBW+1-I 

WRTTF(6.100)X0(TT) ♦UO ( T T ) ♦ WK? ( T ) 




H 

VO 


XD(II)=WK1 (I) 

UD(II)=WK2(I) 

70 CONTINUE 

C*****SM00TH UPSTREAM DATA 
76 IF (N4.EO.O ) 60 TO 99 
99 RETURN 

100 FORMAT ( * XD=*F10.5* UD=*F10.5* UD ( SMOOTHED ) =*F 1 0 .5 ) 

END 

C 

cccccccccccccccccccccceccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


c 

SUBROUTINE LSQSMTH (M.X.Y.NC.C.MS.XS.YS.IC) 

C 

C 

C* LSOSMTH IS A LEAST SQUARES SMOOTHING ROUTINE FOR EITHER A STRAIGHT * 

C* LINE OR A PARABOLIC FIT OF THE INPUT DATA * 

C* * 

C* M NUMBER OF DATA POINTS USED IN THE SMOOTHING * 

C* X * Y 1-D ARRAYS OF COORDINATES OF DATA TO BE SMOOTHED * 

C* NIC NUMBER OF LEAST SQUARES COEFFICIENTS TO BE COMPUTED PLUS 2 * 

C* NC=4 STRAIGHT LINE FIT * 

C* NC=5 PARABOLIC FIT * 

C* C — 1-D ARRAY OF LENGTH N+2 FOR COEFFICIENTS * 

C* MS NUMRER OF POINTS TO BE INTERPOLATED AFTER FITTING DATA * 

C* XS 1-D ARRAY OF LENGTH MS CONTAINING LOCATION OF POINT TO BE * 

C* SMOOTHED * 

C* YS 1-D ARRAY OF LENGTH MS CONTAINING SMOOTHED ORDINATE OF XS * 

C* IC IC=0 COMPUTE COEFFICIENTS * 

C* IC=1 COEFFICIENTS PREVIOUSLY COMPUTED AND SUPPLIED * 

C* 4 

04444444444444444444444444444444444444444444444444444444444444444444444444444444 

c 


DIMENSION X(M) »Y(M) »C(NC) *XS(MS) »YS(MS) 
N=NC-2 

IFCM.LT.NJGO TO 50 
IF(IC.E0.1)G0 TO 21 
SY=SX=SXY=SX2=SX?Y=SX3=SX4=0.0 
DO 3 1 = 1 .M 
SX=SX*X(I) 

SY=SY*Y(I) 

3 CONTINUE 

C<N*1>=SX/M $ C(N*2)=SY/M 

DO 5 1=1. M 
xm=xm-c(N*D 
Y(I)=Y(I)-C(N42) 

5 CONTINUE 
DO 10 1=1. M 
SXY=SXY4X(I)*Y(I) 

SX2=SX2+X(I)*X(I> 

IF f N.EQ.P) GO TO 10 




o o o o 


lo 

o 


SX?Y=SX2Y*X (I)*X (I)*Y(I) 

SX3=SX3+X(I)*X(I)*X(I) 

SX4=SX4*X (I)*X(I)*X(I)*X (I) 

10 CONTI Ml IF 

IF (N.EQ.2) 00 TO 15 

C (3)=(SX?Y*SX2-SXY*SX3) /(SX?*SX4-SX?*SX2*SX2/M-SX3*SX3) 

C <2) = (SXY-C (3) *SX3) /SX? 

C ( 1 ) =-C (3) *SX2/M 
GO TO 1H 
15 C(2)=SXY/SX2 
C(1)=0.0 

19 DO ?0 1=1 »M 
X(I)=X(I)*C(M+1) 

Y ( I ) =Y ( I ) *C (N*2) 

20 CONTI NIIF 

21 IF(MS.EQ.0)GO TO 35 
00 23 1=1 .MS 
XS(I)=XS(I)-C(N*1) 

23 CONTINUE 

DO 30 I = 1 » MS 
YS(I)=0.0 
00 25 J=1 *N 

YS{I)=YS(I)*XS(I ) +C (N+l-J) 

25 CONTI NIIF 
30 CONTINUE 

00 33 1 = 1 * MS 
XS(I)=XS(I)+C(N+1) 

YS(I)=YS(I)+C(N*2) 

33 CONTI NIIF 
35 RETURN 

50 POINT 5 1 f M » N 

51 fodmat (" 1 H ♦/////» * M.LT.N IN LSRSMTH. M=**I3** M=*I3) 

STOP 

END 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

THIS SUBROUTINE TAKES INPUT DATA AMD CALCULATES VELOCITIES AND 
WALL LOCATIONS THEN FILLS APPPOPRIATE ARRAY LOCATIONS. 

SUBROUTINE W ALLOT A (NYPP»YPP» IPRT) 

REAL M.MSET 
INTEGER TEST .RUN. POINT 

COMMON/ JACK/N J *N JT» NJB. X J (100) * YJ ( 100) »UJ ( 1 00) »VJ(10 0) ♦ 

SUJN (100) . VJN (100) . YJN (100) .DYJ(IOO) 

COMMON/P ARAM/HE AD (8) . TE ST » RUN »POI NT ♦ M ,RN » ALPHA » TEMP .KLAX (10)-. 

SMPLAX •PT.PINF »UINF»IAVG»IS(4) • I OPT ( 4 ) «EC(4) .GAMMA «MSF T 
COMMON /CONTOUR /NO . NTW * NDS « MBIfi »NIJS* XD ( 1 0 0 ) » YO ( 1 0 0 ) .UD ( 1 00 ) ♦ VD ( 1 00 ) 

OTMFNSTOM YPP ( MYPP ) • .JOPT ( 2 ) 



C*#* I K 'f E T T OF LIPPFP WALL JACK COORDINATES KOP Y 0 • (VP/UP)* AM;j (VJ/UJ) 

C 

CALI iP \/ 1 ( y J ( 1 ) « y J ( 1 ) * YPP*NJT * 0 • 0 ♦ I OPT ( 1 ) « EC ( 1 ) • * C ( 2 ) • I ) 

On 305 T = 1 • K T v’ 

YP ( I ) = C» W? (xn ( I ) * M JT * X J ( 1 ) fYJ ( 1 ) • YPP*0«0 ) 

VP ( I ) =COPVO ( x o ( I ) ,NJT*XJ(I) ,YJ(1) »YOP»fu0) 

305 CO> T I^'iF 

On 306 T=1*NJT 

VJ(I ) =Ci*RVrUXJ(I) *WJT«*J(1> • Y J ( 1 ) *YPP«0.0) 

3 Oh C P ^T I OOF 
C 

C*****SPLI f *’F ^ T T OF t OWEP WALL JACK COORDINATES POP YOt (VO/UP) * AwO (VJ/UJ) 

r 

CALL HMCi «Pv 1 ( X J (K JT + 1 ) 9 Y J ( N JT ♦ 1 ) • YPP n NJ U « 0 • 0 ? IO^T ( H ) * t C ( 3 ) « *■ C ( A ) « 

^ I r Q D ) 

On 307 I=l.on^ 

I T=K Ti- ♦ vOS* I 

Yf; ( T I ) = n R\/^ ( XD (I I ) t N,* JP * X J { >■ JT ♦ 1 ) ♦ YJ (NJT ♦ 1 ) ♦ Y p P • 0 * i) ) 
vn ( j r ) = c ! pvh ( xo ( J I ) »\>Jp , x j (NJT+ 1 > * Y j ( NJT ♦ i ) * Y°p* o . o ) 

30 7 COr T T r '• iF 

00 3 0’- T = l.^jH 
T T = M J T + T 

m j ( t i ) =ru^V! ( y j ( J I ) 9 NJ^ * X J ( N JT + 1 ) * Y J ( M J T ♦ 1 ) * YP U * o • ! > ) 

3 0> 3 C 0"- ( T I '”jF 
r 

C*tt***C A L CPL ATT HV 0«F (VO/UD) COMPONENT ALCVC THK hO^NST^*:. < ,y *“ K 5 ' 

r 

VOAVOr: ( VO (O'Tv; ) +V0 < NTW+NDQ ♦ 1 ) ) / 2 9 0 
on 3i5 i - 1 • 'v(5 

1 I=Ni T -' + T 

VO (IT) =vna vc 
315 CONTlmjF 

C ... r .. 

C * #C A| _C5 II. AT TOO OF ( VO/UP) CO^PONF NT AL r ‘\*o T Hr. ..VS Tun ' *" 

r 

yPAVO= ( VO(?JT;.+-\*r *VH ( 1) ) /2m 0 

HO 3?o T = l*M.iS 
I T = N T ! n S + ^ *.■» ♦ ] 

VO (IT) =vr A VO 
32 0 COMTTNt!F 

c 

C*tt***C ALCULATTON OF I jf* A NO VP VELOCITIES 

c 

no 3?s i = it r ‘n 

no ( I ) = t * 0 A L ( l ntJP ; F,U') ( I ) fVO ( T ) fOA*- v LA ) 

vn ( I ) =v/o ( I ) *un ( i ) 

3P6 rn^T I \*i »F 


C»****C ALCUL AT TOM OF iij AMD VJ FROM UD DATA 
C 

w JO D T ( 1 ) = 1 

to JOPT(?)=l 

EC1=0.0 
ECM=0.0 
ST = 0 . 0 

CALL HMCMRVl ( X D ( 1 ) * UD ( 1 ) , YPP » MW * ST , JOPT * EC 1 » ECN, t£PP ) 

DO 340 1=1. MJT 

D J ( I ) =CliPV? ( X J ( I ) * NTW , XD ( 1 ) »UD(1) , YPP, ST) 

V J ( I ) =VJ ( I ) *UJ ( T ) 

340 CONTINUE 

IRsMTW+MOS+1 

CALL SWITCH (NRw • XD ( IB ) ) 

CALL SWITCH (N8W»l ID ( IB) ) 

CALL SWITCH (NSW, VD( IB) ) 

JOPT ( 1 ) = ). 

JOPT(?)=l 
FC1 =0 • 0 
FCM=0.0 

CALL HmCHRVI (XD(IR) *UD(IB) , YPP » NRW , ST * JOPT * EC 1 » t CM, IF RP ) 
00 350 1=1, MJB 

UJ (N JT + I ) =Cl)P V2 ( X J (N JT* I ) »NBW*XD( IB) »U0(IB) , YPP, ST) 
v/ J ( M JT+ I ) =V J ( M JT ♦ I ) *U J (MJT ♦ I ) 

350 CONTINUE 

CALL SWITCH (NBW.XO (IB) ) 

CALL SWITCH {MBW,l.iO( IB) ) 

CALL SWITCH (NBW,VD( IB) ) 

C 

C*****REM0V£ FPEESTREAM COMPONENT FROM UD AND UJ TO OBTAIN 
C«****thE PERTURBATION QUANTITIES 
C 

DO 355 1=1, N J 
UJ(I)=UJ(I)-UINF 
355 CONTINUE 

00 360 T = 1 « N D 
UD ( I ) =UD ( I ) -U I NE 
360 CONTINUE 
C 

C*****CL0SF DATA CONTOUR 
C 

ND=NO+ 1 
XD(ND)=XD(1) 

YD (ND) =YD ( 1 ) 
t ID ( ND ) =l!D ( 1 ) 
vd (MD) = vn ( l ) 

c 

C<mhh»*LO0K AT RESULTS OE THIS SUBROUTINE 
C 


o o o 


IFdPPT.FO.DRO TO 390 
WRITE (6*398) l* IMF 
00 365 1=1 ♦ NJ 

WRITE ( 6 • 399) I ♦ X J ( I ) * Y J ( I ) . U J ( I ) ,VJ(I> 

365 CONTIN'JF 

WPITF (6*397) 

00 370 1 = 1 * ND 

WPITF (6*399) I ♦ X T( I ) * YD ( I ) » UO ( I ) * VD ( I ) 

370 CONTIMUF 
390 RETURN 

397 FORMAT { "1 "//* WALL OATA ( SC ALE P ) *//5X* I 

D18X*X0*18X*Y0*1FX*UD*1«X*V0*/) 

398 FORMAT ( "1 "//7X*UI NF=*G1 2» 6//7X* JACK DATA (SCALED) *//5X*T* 
S1PX*XJ*18K*YJ*1«X*UJ*18X*VJ*/) 

399 FORMAT (3X13*4 (5>F15. 9) ) 

END 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccoccccccccccccccccccccccc; 

THIS FUNCTION CALCULATES THE U COMPONENT OF VELOCITY 

FUNCT I ON IlC AL ( M * UI MF * C D ♦ THE T A * GAMMA ) 

REAL M 

T2=l . 0 *THFTA*THFT A 
IF (M.LE.O.O) GO TO 10 
p= (GAMMA-1.0) /GAMMA 
A 1 = 1 . 0- ( 1 ,0 +GAmma*v*M*CP/2»0) **P 

A2=1.0+2.0*Al/( (GAMMA- 1.0 ) *m*M) 

UCAL=UINF*S0PT(A?/T2) 

RFTUPN 

10 A?= 1 . 0— CP 

IF ( A2.LT. 0.0) A2=0.0 
UCAL=UTNF*SORT (A2/T2) 

RFTUPN 
END 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

C THIS FUNCTION REORDERS THE INDICES OF ARRAY X FROM HIGH TO LOW 

c 

SUBROUTINE SWITCH (N*X) 

0 IMF NS I Of! X ( N ) 
no 10 1 = 1 *n 
J=M*1-I 

to IF (J.LE.T) GO TO 11 

w S=X(J) 

X ( J ) = X ( I ) 

X(I)=S 

10 CONTINUE 

11 RETURN 



ooono oooo ooo 


c 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

THIS SUBROUTINE SCALES THE INPUT USING THt GOETHERT SCALING RULES 


SUBROUTINE SCALF (N,X»Y.DY*U.V«m,MSET»IS* IT? IP«T) 

REAL M»»«SET 

OIMFNSIOM X ( N ) »Y (N) *U CM) . V (M) »L'»Y ( N ) 

GETA2=1 .0-M*M 
HETAsSOPT (PETA2) 

RET AMS=SORT ( 1 . 0-NSET*MSFT ) 

*****15=0 COMPRESSIBLE TO INCOMPRESSIBLE 

»e###I5sl INCOMPRESSIBLE TO COMPRESSIBLE 

TF (IS.EO.O. AMD. IPRT.EO.O) WRITE (6*25) 

IF (IS.EO.l. AND. IPRT.EO.O) WRITE (6*26) 

IF (IT.F0.1. AND. IPRT.EO.O) WRITE (F* 2B ) 

IF (IT.EO. 2. AND. IPRT.EO.O) WRITE (6*29) 

IF (IT.EO. 3. AMO. IPRT.EO.O) WRITE (6*30) 

IF (IS.EO.O) GO TO 5 
BETAMS=1 .0/BETAMS 
BET A = 1 . 0/8FTA 
RFTA2=1 .0/RFTA2 
S DO 10 1=1. N 
X(I)=X(I) 

Y (I)=PETA*Y (I) 

TF(IS.EO.l)OY(I)=BETAMS*0Y(I) 

IJ (I ) «BFTA2*U ( I ) 

V ( I ) =9ET A*v ( I ) 

IF ( IT.EO. 1. AMD. IPRT.EO.O) WRITE (6, 27) I * X ( I ) * Y { I ) 
IF(TT.EO.2.AND.IPPT.FO.0) WRITE (fo. 27) I * X ( I ) , Y ( I ) »U ( I ) 

IE (IT.EO. 3. AND. IPRT.EO.O) WRITE (6. 27) I,X(I) ,Y(I) * 0 Y ( I ) *U(I) »V(I) 
10 CONTINUE 
RETURN 

25 FOPMAT(//« COMPRESSIBLE TO INCOMPRESSIBLE SCALING*) 

26 FORMAT <//* INCOMPRESSIBLE TO COMPRESSIBLE SCALING*) 

27 FORMAT (?XI3*5(2X*G12,6) ) 

28 FORMAT (//*X*I*1?X*XJ*12X*YJ*/) 

29 FORMAT (//4X*I*12X*X0*12X*YD*12X*CP»/) 

30 FORMAT ( //4X* I*12X*X J*1 2X*YJ*1 ?X*DY*12X*UJ* 1 2X*V J*/ ) 

END 


THIS SUBROUTINE INTEGRATES THE VELOCITY INTEGRAL TO OBTAIN THE JACK VELOCITY 

SUBROUTINE VJACK(IOPT) 

DIMENSION I OPT ( 4 ) 

r ) 



COPMOW/JACK/UJ.f JT».MJH,XJ(1G0) » Y J ( 1 00 ) .UJ(lUO) .VJ(IOO) , 

‘St I JM (10 0). V JM (100). Y JM ( 1 0 0 ) » D Y J ( ] 0 0 ) 

COmmom/COMTOUR/ND * NTW ♦ MDS * NP w * M IS * XD (100) , YO ( 100 ) .Ml"' (100) . VO < 1 00 ) 
PT=ATAn?(0.0.-1.0) 

T VJ=0 
Nf)Ml=MD-] 
do 40 j=i«mj 

IF ( J.FQ. IVJ) tvRITF (6,96) 

«JM ( J) =0.0 
1 1 JM ( J ) =0 . 0 
00 30 I = J»Mf)Ml 
I5 Tp=TSTRIP ( J. I ) 

IF ( TSTP.PO. J ) CALL STRIP ( J . I . I + 1 ,USTR IP. VSTRIP) 

I c ( ISTR.FQ.2)CALL SI MG ( J, I .USTRIP. VSTRIP) 

VJM ( J) =V JM ( J ) + VSTR I P 
UJM ( J) =IIJM( J) -USTPIP 

TF ( J.EO. IVJ) VRITF (6,99) J, I .XJ ( J) , YJ ( J) ,X0 ( I ) . YO ( I ) . 
*xn(I+l),YD(I+l) . VSTRIP.VJM ( J) * 1STR* UST»IR. UJM ( J) 

30 COMTImuF 
°IDIV=dI 

IP (IPORY(J) .EO.?)PIDIV=2.0*PT 
VJM ( J) = VJM ( J) /Pimv 
UJM(J) = UJN( J) /PTPIV 
4 0 COMTIMI.IF 

IF(I0PT(1) .EO.llfiO TO 42 
VJN(1)=VJ(1) 

IIJM(l)=tlJ(l) 

42 IF(IOPT (?) .FO.l )GO TO 44 
VJM (UJT) =VJ (NJT) 

UJN (NJT) =IIJ (NJT) 

44 I F ( I OPT (3). FO.l) GO TO 46 
VJM (MJT* 1 ) =V J (MJT+ 1 ) 

UJM (M JT* 1 ) =UJ (M JT* 1 ) 

46 TF(TOPT(4) .F.O.l ) GO TO 48 
VJM (MJT+MJR) =VJ (MJT+NJB) 

UJM (NJT+NJB) =UJ (OJT+NJB) 

48 WRITF (6.97) 

00 50 1=1, MJ 
OVJ=V JM ( I ) — V J ( I ) 

OUJ=U JM ( I ) -U J ( I ) 

■VPITF(6.Rfl) I.XJ(I) ,YJ(I) »UJ(I ) , U JM ( I ) ,VJ(I) ,VJM(I) .OVJ.I'UJ 
50 CONTINUE 
RETURN 

96 FOPMAT (" 1") 

97 FORMAT Cl”//* MFW JACK PPJRTURHAT IOM VELOCITY CALCULATIONS* 

S* (SCALFO)*//* I* 

< R6X*XJ*3X*YJ*9X*UJ*6X*UJN*5X*VJ*RX*VJN*9X*0VJ*9X*UUJ*/) 

98 FORMAT ( 13, 2 (?XF 8.4) .2 ( 2XF7 . 2 ) , 4 ( ?XG 1 0 .4) ) 

99 Format (* J=*I2* I=* I 2* XJ=*F7.3* YJ=*F8.4* XO=*F7.3* Y0=*F8.4 
5* X0=*F7.3* YO=*F 8.4* VSTRIP=*G 12.6* VJN**S12.6* ISTPIP=*I2/ 

f 7 9 X * I(STPTP = *G1 P.6* |(JM = *G) 2.6) 



o o o o 


END 


c 

cccccccccccccrcccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

THXS Sll«POl)TINF CALCULATES THF JACK VELOCITY COMPONENT RESULT INC- FROM 
sr»MF ST»T P EXCLUDING THAT CONTAINING THE SINGULARITY 


SUBROUTINE STRIP ( I • K * L » USTR I» * VSTW IP ) 

COMMON/ JACK/NJ»NJT» NJR, X J ( 1 no ) * YJ ( 100 ) ,UJ ( 100 ) , V J ( 1 00 ) , 

■SljJN ( 100 ) . VJN (100) »YJN ( 100 ) ♦ DY J ( 1 0 0 ) 

COMmON/CONTOUR/MD ,NTW * NOS , NR* » NUS » XD ( 100 ) » YD ( 100 ) *UD (100) , VD (100) 
XL = XD (L ) -XD ( K ) f. YL=YO(L)-YD(K) S SL=SOPT ( XL*XL + YL*YL ) 
RX=XJ(I)-XD(K) « RY=YJ(I)-YO(K> 

UOIF=UD (L) -UP (K) S VDIF=VD (L ) -VD (K) 

Al=-9X*llD(K) % A?=-PX*VO (K ) f A3=-RY*UO(K) * A4=-BY*VD(K) 
ni = (-9X<NIDTF*XL*U0 (K) ) /SL 6 B?= (-RX*vniF+XL*VD (*) ) /SL 

R3= (-RY*IIDIF+YL*UD (K ) ) /SL y R4= ( -RY*VPI E+ YL*V0 (* ) ) /SL 
SL2=Sl*SL 

Cl=Xl AU0TF/SL2 " C2=XL*VDIF/SL? 

C3=YL*UDIF/SL? <S C4=YL*VDTE/SL2 

01=RX»RX*RY*RY S 02= -2* (RX*XL+RY*YL) /SL 

F 1 = XL /St. S E2=YL/SL 

A0= A 1 -A4 S RD=B 1-B4 * CP=C1-C4 

AP=A3*A2 S RP=B3+R2 * CP=C3+C2 

0=4. *01-02*0? 

IF(0)630»610,605 
60S 0RT=S0RT(0) 

T1 =?• * AT AN ( SL *QRT / ( 2 . 0*01 +SL*D2 ) ) /OPT 
APG=(01+02*SL+SL?) /01 
I F ( ARG . L F . 0 • 0 ) GO TO 650 
T?=.5*AL0G ( APG) -.5*0?*T1 
T3=SL~0?*T?-01*T 1 
GO TO 620 

610 Tl-2. 0*SL/ ( 02* ( 02/2. 0*SL ) ) 

ARG=2.* (02/2. 0*SL) /02 
IF (ARB.l.F.O.O)GO to 650 
T?=AL0G(ARG)-.5*02*T1 
T3=SL ”0?*T2- (02/2. ) * (02/?. >*T1 
620 VSTRIP= (FI* ( AD*T1 +R0*T?+C0*T3 ) +E?*(AP*T1+BP*T2+Cp*T3) ) 

USTRIP= (El* ( AP*T1+RP*T2*CP*T3) -F2* ( AD*T1+B0*T?+C0*T3) ) 

RETURN 

630 PRINT 640*0 . 

640 FORMAT (//,*FRROP- — THE VALUE OF 0 IS NEGATIVE, Q-* * F 1 5 . « * // ) 

STOP 

650 PRINT 651 » I »K *L * APG 

651 FORMAT!/,* I=*»I2»* J=*»I?»* J+1=*,I?»* ARG=*»F?5, 15, 

*/,* F»poR— -ARGUMENT OE LOG IN STRIP. LT.O*) 

STOP 

FNO 



to 

-J 


cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c THIS SUBROUTINE EVALUATES THE VELOCITY COMPONENT ALONG THE SINGULAR STRIP 


SUBROUTINE SING(I,J,USING*VSING) 

COMMON/JACK/NJ»NJT»NJB t XJ(100) *YJ(100) ,UJ(100) ,VJ(100) * 

$UJN (100) * VJN (10 0) ♦YJN(IOO) ♦OYJ(IOO) 

PC-SO^M^?T? / ^^y^*o S ’ NBW,NUS,XD(100) »YO(100),UD(100),VO(100) 
PC-SORT ( (XJ(I)-XO(J) ) # *H+(YJ(I)-YD(J) )**?) 

XJX1=XJ(I)-XD(J) 

YJY1SXJX1*(YD(J*1)-YD(J) )/(XO(J*l)-XD(J) ) 

PC=SORT ( YJY1*YJY1+XJX1*XJX1 ) 

SL = SQRT ( (XD(J*1)-XD(J) )**2*(Y0(J+1)-Y0(J) )*«•?) 

APG=SL/PC-1 . 0 
IF (ARG.LE ,0« ) GO TO 750 
T1=UD ( J+l ) -U0 ( J) 

T2=UD (J)+T1*PC/SL 
T3=AL0G(ARG) 

T4=VD(J*1)-VD(J) 

T5=VD ( J) +T4*PC/SL 
VSING=T1+T2*T3 
USING=T4-fT5*T3 
RETURN 

750 PRINT 751 * I * J* ARG 

751 FORMAT (* I=*»I2»» J=**I2»* 

SARGUMENT OF LOG IN SING LE 

STOP 
END 


ARG=*,F25.15,/»* ERROR — 
0») 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


c 

C 

c 


Jow«[ U op T s?ngSur R s"ip " hether ™ e strip seing “nsidered is A 


c 


FUNCTION ISTRIP ( I J» ID) 

COMMON/ J A CK/NJ * N JT ♦ NJB ♦ X J ( 1 00 ) » Y J ( 1 0 0 ) ♦ U J ( 1 00 ) 
SUJN(IOO) ♦ VJN (100) » YJN (100) »DYJ(100) 
COMMON/CONTOUR/ND *NTW»NDS»N8W »NUS »XD(100) *YD(1 


»VJ (100) * 

00) » UD (10 0) fVD(lOO) 


C*»***ISTRIP=1 NONSINGULAR STRIP 

C**** # ISTRIP=2 SINGULAR STRIP 

C 


I STRI P= 1 
SJ=XJ(IJ) 
TJ=YJ(IJ) 
S1=XD (ID) 
T1=YD ( ID) 
S2=XD (ID+1 ) 
T2=YD ( ID+1 ) 
niF=S?-Sl 



DODO 


tsj 

00 


PR0D=T1*TJ»T? 

10 {Jl"jle?!M’.wn.SJ.LT.Sl).»Mn. (PPOO.LT. O.ODISTBIP-2 

20 IF USJ.OT.Sl . AND.SJ.LT.S2) .AMO. (PROD. (IT. 0 . 0 ) ) ISTRIP=2 
30 CONTINUE 
RETURN 
END 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

THIS FUNCTION DFTERMINES WHETHER A POINT IS A POINT ON THE DATA SURFACE 
OR A FIELD POINT 

COMMON/ JACK /NJ ♦ N JT . NJB » X J ( 1 00 ) » Y J ( 1 00 ) ♦ UJ < 1 00 ) ♦ V J ( 1 00 ) * 

Q 

C*****IP0RY=1 - ( XJ ♦ Y J ) ON OATA SURFACE 
C*****ibDRY=2- — ( X J * Y J ) IS FIELD POINT 

C 

IBDRY=1 

IFUSjlLT.XO(l) ) .OR. (SJ.RT.XD(NTW) ) ) IBDRY=2 

RETURN 

END 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c THIS SUBROUTINE INTEGRATFS THE DELTA THETA VALUES TO OBTAIN NEW 

c WALL LOCATION PREDICTIONS 

C SUBROUTINE DELYJ(NTHETA*THETA,NYPP*VPP.UINF) 

htmfmqtom THETA (NTHETA.) t YPP (NYPP) * JOPT (2) 

COMMON /J ACK/NJ * NJT .NJB* X J (100) *YJ(100) .UJ ( 100) f VJ ( 100) * 

$UJN(100) .VJN(IOO) .YJN(IOO) tDYJ(lOO) 

JOPT ( 1 ) =1 

JOPT (2) =1 

EC1=0 . 0 

ECNsO.O 

DO 10 1=1 »NJ 

||JN ( I ) =UJN ( I ) ♦UINF 

UJ ( I ) =UJ ( I ) ♦UINF 

THETA ( I ) =V JN ( I ) /U JN ( I ) 

YJN(I)=VJ(I)/UJ(I) 

10 CALL^HMCURVl (XJ. THETA. YPP.NJT. 0.0. J0PT.EC1.ECN.IERR) 

no IS T = 1 .NJT 



no non non 


r>YJ(I)=CHPVI (XJ() ) *XJ(I) ,NJT«XJ,THETA,YPp,0.0) 

IS CONTINUE 

CALI HMCI.IPV1 (XJ,YJN,YPP«NJT»0.0,JOPT*EC1 ttCN. IERR) 

00 ?0 1=1, NJT 

OYJ ( I ) =f>YJ ( I ) -CUPVI (XJ(1) , X J ( T ) ,NJT,XJ,YJN,YPP,0.0) 

?0 CONTINUE 
N1=NJT + 1 
N2=NJT*HJR 

CALL H^C'PVl (XJ(N1) , THET A (Ml) , YPP ♦ N JE , 0 . 0 , JOPT , EC 1 , PCM » 1 F RR ) 
no 25 l=wi,N2 

OYJ(I)=ruRVl (XJ(N1) »XJ(I> ,NJH,XJ(N1) ,THET4(M1) , YPP, 0.0) 

25 CONTINUE 

CALL HMCMRVl (XJ(N1) *YJN(N1) , YPP , NJB , 0 . 0, JOPT , EC1 , ECN , IEPR ) 

00 30 I = N 1 » N 2 

OYJ(I)=r»YJ(I)-CllRVI ( X J ( N 1 ) ,XJ(I) , NJb , X J ( N 1 ) ,YJN(N1) , YPP, 0.0) 

30 CONTINUE 
WRITE (6,08) 

00 40 1=1, NJ 

WRITE (6,09) I, XJ(I) ,YJ(I) tUJ(T> ,UJN ( I ) ,V/J ( I ) , V JN ( I > , Y JN ( I ) , ThETA ( I ) 

*,DYJ( I ) 

40 CONTINUE 

RMSDYT=RMS(NJT»DYJ) 

PMSf)YB=PNS(NJB,OYJ (NJT+1 ) ) 

WPITF (6,97) RMSDYT.RMSDYR 
RETURN 

97 FORHAT(//* RMS(OYJ(TOP) )= *01?. 6* 3x5 ( OY J ( ROT ) ) = *612.6) 

96 FORMAT ("1 M //* PREDICTED WALL CORRECTION (SCALED)*// 

, S2X*I*6X*XJ*8X*YJ*9X*UJ*7X*UJ\ | *5X*VJ*9X*VJN*9X*THETA*6X 
$*THET AN*tvX*DY J*/) 

99 FORMAT (T3,?(?XE8.4) , 2 ( 2XF7 . 2 ) , S ( ?XG 1 0 . 4 ) ) 

END 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

THTS FUNCTION STATEMENTS DETERMINES THE RMS VALUE OE ARRAY A 

FUNCTION RMS (N, A ) 

DIMENSION A (N) 

RMS=0.0 
DO 10 1=1 ,N 
PMS=RMS*A (I)*A(I) 

10 CONTINUE 

rmS=SORT (RMS/N) 

RETURN 
END 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


SUBROUTINE HMCURVl ( X , Y , YPP , N , ST , OPT , EC 1 »ECN» IEPR) 


THTS CUHROUTTNF FTTS A ClIRTC OR TENSION SP) TNF THROUGH A SET 



o o o o o r> 


OF Y VS X VALUES, 


AUTHOR — HARRY MORGAN NAS A /L ANGLE Y LS AH/ST AD 

(MODIFIED BY JOEL EVERHART NAS A/L ANGLE Y ) 


OFSCR I RT I n?> J OF INPUT AND OUTPUT FOR SUBROUTINE SPLINE 

DESCRIPTION 


C* 

C* 

C* 

C* PAPAMFT FR 
C* 

Y - INPUT ARRAY OF I NOFPF N’DE NT VARIABLE. 

INPUT ARRAY OF oepfndfnt variable* 
n U M v Y WORK ARRAY OF DIPENSION 10*M ON INPUT. 

ON OUTPUT * ARRAY OF SECOND DERIVATIVES. 

NUMBER OF INPUT values OF X AND Y. 
spline TENSION (ST=0.0 same as CUBIC SPLINE). 

END CONDITION OPTION ARRAY OF LENGTH 2. OPT (I) 
CORRESPONDS TO CONDITION AT X(l)* OPT (2) CORRESPONDS 
TO CONDITION AT X <N) . 

1 - FC1 OR ECN UNKNOWN 
? - FCl-YP(l) OP ECN=YP(N) 

3 - F C 1 = Y PP ( 1 ) OR ECN=YPP(N) 

END CONDITION AT END OF Y VS X CURVE. 

ERROR CODE 
0 - NO ERROR 
1 - N IS LESS THAN 3 
? - X I S NOT MONOTONICALLY INCREASING 


C* 

C* 

C* 

C* 

C* 

c# 

c* 

C* 

c* 

c* 

c* 

c* 

c* 

c# 

c* 

c* 

c* 


Y 

Y pp 

N 

ST 

OPT 


FCN 

IFPR 


C* 

C* RESTRICTIONS 

C* X MUST R F NONOTON I C ALL Y INCREASING 

C* M MUST be greater THAN OR EQUAL TO 3. 

C* X AND Y MUST BE DIMENSIONED BY AT LEAST N T K CALLING 

C* PROGRAM. 

0* ST SHOULD BE NO GREATER THAN SO. 

C* 

C************************************* **************** ****************** 

c 

DIMENSION X(N)f Y ( N ) * YPP ( 1 ) • 0PT(2> 


* 

* 

* 

* 

# 

■it 

* 

* 

* 

■it 

■* 

•it 

* 

•* 

* 

«t 

# 

■it 

■it 

•it 

*t 

■it 

•it 

•it 

■it 

•it 

«* 


C 

C 


C 

c 

C 


INTEGER OPT 

STNH (X) =0.5* (EXP (X) -EXP <-X) ) 

COSH (X) =0.5* (FXP(X) ♦FXP(-X) ) 

CHECK ON MONOTONI CALL Y INCREASING X 

lFpp=n 

TE (N.l T.3) GO T 0 3? 



ooo ooo o o o o o o o o o 


00 1 I = ? . M 

IF (X(I).LF.X(I-l) ) GO TO 33 
1 CPUT J MI.IP 

C OMPHTF Tt MS I ON PARAMETER T 


T=sST*FLOAT (\‘-l )/(X(N)-X(l)) 


COMPUTE NEEDFH PROGRAM IMFGFP VALUES 

m=n-1 

Ml=M+l 
U? = K<1 + M 

m3=M2+n 
M4 = U 3 + U 
N 5 = N 4 + "-J 
N6 = M5*N 

r7 = N6 + i\i 
M8 = M7 + m 
M 9 = K' 8 + M 

LOAD MATRIX ROWS 2 THRU U-l 

00 ? 1=1 tM 

H2 = X ( I + 1 ) -X ( I ) 

IF (ST. IF. 0.0) YPP(I)=H2/6. 

IF (ST.GT.0.0) YPP ( I ) = ( 1 ./H2-T/SINH (T*H?) )/(T*T) 

K=2*N+ I 

YPP ( K ) =YPP ( I ) 

IF (I.EO.l) GO TO 2 
K = K + N 

IF (ST.LF.0.0) YPP (K)=(H?+Hl)/3# 

IF (ST.GT. 0 . 0 ) YPP (K) = < T*COSH (T*H1 > /SINH < T*H1 > -1 ./Hi *T*COSH ( T*H? ) /SI 
SNJH(T*H?)-1./H2) /(T*T) 

K = K + N 

YPP(K)=(Y(I+1)-Y(I) ) /H2- ( Y ( I ) -Y ( 1-1 ) ) /Hi 

2 H 1 = H? 

LOAD matrix ROWS 1 AMO N ( FMD CONDITIONS) 

YPP (M3) =YPP (N4-1 ) =1 .0 
YPP { N? ) =YPP (M) =0 .0 

POW 1 

GO TO (3*4*6) 0 D T (1) 

3 A 1 =X (?) -v (3) 

A 2= X ( 3 ) -X ( I ) 

A3 = X ( 1 ) -X (?) 

YPP ( N4 ) =? • * ( Y ( 1 ) * A I + Y ( 2 ) * A2 + Y ( 3 ) * A3 ) / ( A I * X ( 1 ) A?* x ( 2 ) **2+ A3* 

SX (3) **?) 



o o o o o o 


GO TO 7 

4 H1=X(2)-X(1) 

YPP (N4)=EC1— (Y(2)-Y(l) >/Hl 
IF (ST.GT.0.0) GO TO 5 
YPP (N3) =-Hl/3. 

YPP (M?) =-Hl/6. 

GO TO 7 

5 YPP ( N3) =<1./(H1*T) -COSH (T*H1)/SINH(T*H1) ) /T 
YPP (N2) = ( 1 ,/SINH (T*H1 ) -1./ (H1*T) ) /T 

GO TO 7 

6 YPP ( N4 ) =FC1 

ROW N 

7 GO TO ( ft *9* 1 1 ) OPT (2) 

8 A1=X (N-l ) -X (N) 

A2=X ( N ) -X (N-2) 

A3=X(M-?)-X (N-l ) 

YPP (N5-1 ) =2.*(Y(N-2) *A1+Y (N-l ) #A2*Y (N) *A3) / ( A1*X (N-2) **2*A2*X (N— 1 ) 
*»*2+A3*X (N) **2) 

GO TO 1? 

9 H1=X (N) -X (M) 

YPP (NS- 1 )=ECN-(Y(N)-Y(M) )/Hl 
IF (ST.GT.0.0) GO TO 10 
YPP (M) =Hl/6. 

YPP (N4-l)=Hl/3. 

GO TO 1? 

10 YPP ( M ) = ( 1 ./ ( H 1 *T ) -1 ,/SINH ( T*H 1 ) ) /T 

YPP (N4-1 ) = (COSH (H1*T) /SINH (T*H1 ) - 1 . / ( H1*T ) ) /T 
GO TO 1? 

11 YPP (N5-1 ) =ECN 

SOLVE TRIDIAGONAL MATRIX 

12 YPP (N7) =YPP (N3) 

YPP (N8) =YPP (N4) /YPP (N3) 

YPP (N9) =YPP (N2) /YPP (N3) 

DO 13 I=2»N 
11=1-1 
12=1-2 

YPP (N7+I1)=YPP(M3+I1)-YPP(I1)*YPP(N9*I2) 

IF (I.EO.N) GO TO 13 
YPP(N9+I1)=YPP(N2+I1)/YPP(N7+I1) 

13 YPP ( N8+ 1 1 ) = (YPP (M4+I1 ) -YPP (ID #YPP(N8+I2) ) /YPP (N7 + 1 1 ) 

YPP (N6-1 ) =YPP (N9-1 ) 

DO 14 1 = 1, M 
1 1 =N6-1 -I 
I 2=N- I 
T3=T2-1 



noooooooooooooaoooooo 


\ 


c 

c 

c 


U 


15 


YPP (I1)=YPP(N8+I3)-YPP(N9+I3)*YPP(I1*1) 
K = 0 

00 15 1=1, N 
11 = 1-1 

YPP (N6+K ) = YPP (N5+ 1 1 ) 


K=K ♦ 1 

YPP(M)=YPP(N3-2) 


REARRANGE SOLUTION IN YPP ARRAY 


c 


c 


00 16 1=1, N 
1 1 =N6 ♦ I - 1 
16 YPP ( I ) =YPP (II) 


RETURN 

32 IERR=1 
RETURN 

33 IFRR=2 

WR I TF ( 6 , 35 ) (I ,X ( I ) »Y(I) » 1=1 ,N) 

RETURN 

35 FORMAT (1H1,5X*33HX IS NOT MONOTONICALLY 
$X,19X,1MY/(5X,I5,5X,E15.6,5X,F15.6) ) 

END 


INCREASING//9X»1HI»13X,1H 


cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


FUNCTION CUPV2 (T»N,X,Y,YP, SIGMA) 
INTEGER N 

REAL T,X (N) ,Y (N) ,YP <N) , SIGMA 


FROM THE SPLINE UNDER TENSION PACKAGE 
CODED BY A. K, CLINE AND R. J. RENKA 
DEPARTMENT OF COMPUTER SCIENCES 
UNIVERSITY OF TEXAS AT AUSTIN 

THIS FUNCTION INTERPOLATES A CURVE AT A GIVEN POINT 
USING A SPLINE UNDER TENSION, THE SUBROUTINE CURVl SHOULD 
BE CALLED EARLIER TO DETERMINE CERTAIN NECESSARY 
PARAMETERS, 


ON INPUT — 


T CONTAINS A REAL VALUE TO BE MAPPED ONTO THE INTERPO- 
LATING CURVE. 

N CONTAINS THE NUMBER OF POINTS WHICH WERE SPECIFIED TO 
DETERMINE THE CURVE. 

X AND Y ARE ARRAYS CONTAINING THE ABSCISSAE AND 
ORDINATES. RESPECTTVFI Y, OF THF SPFCTFTFD POINTS. 



non ooo ooooooonoooooooooooooooo 


YP IS AN ARRAY OF SECOND DERIVATIVE VALUES OF THE CURVE 
AT THE NODES. 

AND 

SIGMA CONTAINS THE TENSION FACTOR (ITS SIGN IS IGNORED). 

THE PARAMETERS N* X. Y» YP» AND SIGMA SHOULD 3E INPUT 
UNALTFRFD FROM THE OUTPUT OF CURV.l . 

ON OUTPUT — 

CURVE! CONTAINS THF INTERPOLATED VALUE. 

NONE OF THE INPUT PARAMETERS ARE ALTERED. 

THIS FUNCTION REFERENCES PACKAGE MODULES TNTRVL AND 
SNHCSH. 


DETERMINE INTERVAL 

I Ml = INTRVL (T.X.N) 
I = IM1+1 


OENORMALIZE TENSION FACTOR 

SIGMAP = ABS(SI6MA)*FL0AT ( N— 1 ) / ( X (N)-X(l) ) 


SET UP AND PERFORM INTERPOLATION 


DELI = T-X(IMl) 

0EL2 = X(I)-T 
DELS = X<t)-X(IM1) 

SUM = (Y(I)*DEL1*Y(IM1)*DEL2)/DELS 
IF (SIGMAP ,NE. 0.) GO TO 1 

CURV2 = SUM-DEL1*DEL2*(YP(I)*(DEL1+DELS)+YP<IM1)* 
* (DEL2+DELS) )/(6.*DELS) 


RETURN 

DEL PI = SIGMAP*(DELl*DELS)/2. 

DELP2 = SIGMAP* (DEL2+DELS)/2. 

CALL SNHCSH (SINHM1 ,DUMMY*SIGMAP*DEL1 »-l ) 
CALL SNHCSH ( S I NHM2 » DUMMY » S IGMAP*DEL2 ♦ -1 ) 
CALL SNHCSH (SINHMS»DUMMY*SIGMAP*DELS»-i ) 
CALL SNHCSH (SINHP1 ♦DUMMY»SIGMAP*DELl/2. ,-l ) 
CALL SNHCSH ( S I NHP2 * DUMMY » S IGMAP*DEL2/2 . » - 1 ) 
CALL SNHCSH ( DUMMY . COSHP1 , DELP 1 * 1 ) 



CALL SNHCSH (DUMMY, C0SHP?*DELP2* 1 ) 

CURV? = SUM+(YP(I)#(SINHM1*DEL2-DEL1#(2.MC0SHP1 + 1.)# 

* SINHP2+SIGMAP*C0SHP1*0EL2> ) + YP ( I Ml ) * (SINHM2* 

* nELl-0EL2*(2.*(C0SHP2+l. ) *SINHP1*SIGMAP* 

* C0SHP2*DEL 1 ) ) )/<SIGMAP*SI6MAP*DELS*(SINHMS+ 

* SIGMAP*OELS) ) 

RETURN 

END 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCGCCCCCCCCCCCCCCCCCCCCCCC 

c 

FUNCTION CURVD ( T * N * X * Y » YP » S I GMA ) 

c 

INTEGER N 

REAL T*X(N) *Y(N) ,YP(N> , SIGMA 
C 

C FROM THE SPLINE UNDER TENSION PACKAGE 

C CODED BY A. K. CLINE AND R. J. RENKA 

C DEPARTMENT OF COMPUTER SCIENCES 

C UNIVERSITY OF TEXAS AT AUSTIN 

C 

C THIS FUNCTION OIFFERFNTIATES A CURVE AT A GIVEN POINT 
C USING A SPLINE UNDER TENSION. THE SUBROUTINE CURVl SHOULD 
C BE CALLED EARLIER TO DETERMINE CERTAIN NECESSARY 
C PARAMETERS. 

C 

C ON INPUT — 

C 

C T CONTAINS A REAL VALUE AT WHICH THE DERIVATIVE IS TO BE 
C DETERMINED. 

c 

C N CONTAINS THE NUMBER OF POINTS WHICH WERE SPECIFIED TO 
C DETERMINE THE CURVE. 

c 

C X AND Y ARE ARRAYS CONTAINING THE ABSCISSAE AND 
C ORDINATES* RESPECTIVELY* OF THE SPECIFIED POINTS. 

C 

C YP IS AN ARRAY OF SECOND DERIVATIVE VALUES OF THE CURVE 
C AT THE NODES. 

C 

C AND 
C 

C SIGMA CONTAINS THE TENSION FACTOR (ITS SIGN IS IGNORED). 

C 

C THE PARAMETERS N» X* Y* YP* AND SIGMA SHOULD BE INPUT 
C UNALTERED FROM THE OUTPUT OF CURVl. 

c 

C ON OUTPUT— 

C 

C CURVD CONTAINS THE DERIVATIVE VALUE. 

C 



OOOOOOO 0-000 ooo ooo ooooooooo 


NONE OF THE INPUT PARAMETERS ARE ALTERED. 

THIS FUNCTION REFERENCES PACKAGE MODULES INTRVL AND 
SNHCSH. 


DETERMINE INTERVAL 

IM1 = INTRVL(T.X.N) 

I = I M 1 + 1 

DENORMALIZE TENSION FACTOR 

SIGMAR = ARS(SIGMA)*FL0AT(N-1)/(X(N)-X(1)) 

SET UP AND PFRFORM DIFFERENTIATION 

DFL1 = T-X(IMl) 

DFL2 = X ( I ) “T 
DELS = X(I)-X(IM1) 

SUM = (Y(I)-Y(IMl ) J/DELS 
IF (SIGMAP .NE. 0.) GO TO 1 

CURVD = SUM+(YP(I)*(?.*DEL1*DEL1-DEL?*(DEL1+DELS) )- 

* YP(IM1)*(2.*DEL2*0EL2-DEL1*(DEL2*DELS) ) ) 

* / (6.*0ELS) 

RETURN 

1 CALL SNHCSH ( DUMMY » C0SHM1 »SIGMAP*DEL1 * 1 ) 

CALL SNHCSH (DUMMY»C0SHM2,SIGMAP*DEL2* 1 ) 

CALL SNHCSH ( SINHMS » DUMMY » SI GMAP*DELS »-l ) 

CURVD = SUM* (YP ( I )* (DEL S*SIGM AP*C0SHM1— SINHMS) ~ 

* YP(IM1)*(DELS*SIGMAP*C0SHM2-SINHMS) )/ 

* <SIGMAP«SIGMAP*DELS*(SINHMS+SIGMAP*DELS) ) 
RETURN 

END 


FUNCTION CURVI <XL»XU*N*X»Y»YP*SIGMA) 

INTEGER N 

PEAL X L * X U * X (N) »Y (N) ,YP (N) » SIGMA 

FROM THE SPLINE UNDER TENSION PACKAGE 
CODED BY A. K. CLINE AND R. J. RENKA 
DEPARTMENT OF COMPUTER SCIENCES 
UNIVERSITY OF TEXAS AT AUSTIN 

THTS FUNCTION INTEGRATES A CURVE SPECIFIED BY A SPIIME 



ooo non oooooooooooooooooooooaoooooooooooooooooo 


UNDER TENSION BETWEEN TWO GIVEN LIMITS. THE SUBROUTINE 
CURV1 SHOULD BE CALLED EARLIER TO DETERMINE NECESSARY 
PARAMETERS. 

ON INPUT — 

XL AND XU CONTAIN THE UPPER AND LOWER LIMITS OF INTE- 
GRATION, RESPECTIVELY. (SL NEED NOT BE LESS THAN OR 
EQUAL TO XU, CURVI (XL, XU,...) .EG. -CURVI (XU.XL,...) >. 

N CONTAINS THE NUMBER OF POINTS WHICH WERE SPECIFIED TO 
DETERMINE THE CURVE. 

X AND Y ARE ARRAYS CONTAINING THE ABSCISSAE AND 
ORDINATES, RESPECTIVELY, OF THE SPECIFIED POINTS. 

YP IS AN ARRAY FROM SUBROUTINE CURVI CONTAINING 
AT THE NODES. 

AND 

SIGMA CONTAINS THF TENSION FACTOR (ITS SIGN IS IGNORED). 

THE PARAMETERS N» X, Y, YP, AND SIGMA SHOULD BE INPUT 
UNALTERED FROM THE OUTPUT OF CURVI. 

ON OUTPUT — 

CURVI CONTAINS THE INTEGRAL VALUE. 

NONE OF THE INPUT PARAMETERS ARE ALTERED, 

THIS FUNCTION REFERENCES PACKAGE MODULES INTRVL AND 
SNHCSH. 


STATEMENT FUNCTION FOR COEFFICIENT ASSOCIATED WITH 
DERIVATIVE TERMS 

TERM (CMM1 ,CMM2,T) = (CMM1-CMM2-SIGMAP*T*SS) / (SIGMAP* 
* SIGMAP*SIGMAP*(SS*SIGMAP*DELS) ) 

DENORMALIZE TENSION FACTOR 

SIGMAP = ABS(SIGMA)*FLOAT(N-l)/(X(N)-X(l) ) 

DETERMINE ACTUAL UPPER AND LOWER BOUNDS 

XXL = XL 
xxu = XII 



non oooooo 


u> 

00 


SSIGN = 1. 

IF (XL .IT. XU) GO TO 1 
XXL = XU 
XXU = XL 
SSIGN = -1. 

IF (XL .GT. XIJ) GO TO 1 

RETURN 7ER0 IF XL .EO. XU 

CUR VI = 0. 

RETURN 

SEARCH FOR PROPER INTERVALS 

1 ILM1 = INTRVL (XXL.X.N) 

IL = ILM1+1 

IUM1 = INTRVL (XXU.X.N) 

IU = IUM1+1 

IF (IL .EQ. IU) GO TO 8 
INTEGRATE FROM XXL TO X(IL) 

SUM = 0. 

IF (XXL .EQ. X ( IL ) ) GO TO 3 
DELI = XXL-X(ILMl) 

0EL2 = X ( IL ) “XXL 
DELS = X(IL)-X(ILM1) 

T 1 = (0EL1+0ELS)*DEL2/(2.*DELS) 

T2 = DEL2*DEL2/ (2.*DELS) 

SUM = T 1 *Y (IL) +T2*Y (ILM1) 

IF (SIGMA .EQ. 0.) GO TO 2 

CALL SNHCSH ( DUMMY ♦ C 1 . S IGMAP*DEL 1 » 2 ) 

CALL SNHCSH (DUMMY»C2.SIGMAP*DEL2»2) 

CALL SNHCSH (SS.CS.SIGMAP#DELS»3) 

SUM = SUM*TERM(CS,C1.T1)*YP(IL> 

* +TERM(C2»0.»T2)*YP(ILM1) 

GO TO 3 

2 SUM = SUM-T1*T1*0ELS*YP(IL) /6. 

* -T2* (DEL 1* (DEL2+DELS ) ♦DELS*DELS) *YP(ILM1)/12. 

INTEGRATE OVEP INTERIOR INTERVALS 

3 IF (IU-IL .EO. 1) GO TO 6 
ILP1 = IL+ 1 

DO 5 I = ILP1 * IUM1 
DELS = X(I)-X(I-1> 

SUM = SUM-*- (Y (I ) +Y (1-1) ) *DELS/2. 

IF (SIGMA .EQ. 0.) GO TO 4 

CAI I SNHCSH (SS.CS.SIGMAP«OFI S.3) 



non 


SUM = SUM*(YP(I)+YP(I-1) )*(CS-SS*SlGMAP*0ELS/2.)/ 

* (SIGMAP*SIGMAP*SIGMAP* (SS+SIGMAP»DELS) ) 

GO TO 5 

4 SUM = SUM- ( YP ( I ) +YP ( 1-1 ) ) *DELS*DELS*DELS/24 • 

5 CONTINUE 
C 

C INTEGRATE FROM X(IU-l) TO XXU 

6 IF (XXU .EO. X(IUMl)) GO TO 10 
0EL1 = XXU-X(IUMl) 

DEL? = X(IU)-XXU 
DELS = X(IU)-X(IUM1) 

T 1 = 0EL1*DEL1/ (?.*DELS) 

T2 = (DFL?+DELS)*DEL1/(2.*0ELS) 

SUM = SUM+T1*Y ( III) + T2*Y < IUM1 ) 

IF (SIGMA .EG. 0.) GO TO 7 
CALL SNHCSH (DUMMY » Cl »SIGMAP*DEL 1 « 2 ) 

CALL SNHCSH ( DUMMY *C2 *SIGMAP*DEL2 * 2 ) 

CALL SNHCSH ( SS ♦ CS »SIGMAP*DELS * 3 ) 

SUM = SUM ♦ TERM ( Cl *0»*T1) # YP(II!) 

* +TERM(CS»C2*T2)*YP(IUM1) 

GO TO 10 

7 SUM = SUM-Tl* (OFl.2* ( DE L1+ DELS) + DELS* DELS) *YP(IU)/12. 

* -T2*T2*DELS*YP(IUMl)/6. 

GO TO 10 

INTEGRATE FROM XXL TO XXU 

8 DELHI = XXU-X(IUMl) 

DELU2 = X( IU> -XXU 
DELL1 = XXL-X ( I U M 1 ) 

DELL2 = X ( I U ) — X X L 
DELS = X(IU)-X(IUM1) 

DELI = X XU— X XL 

T 1 = <DFLU1+DELL1)*DELI/<2.*DELS) 

T2 = (DFLU2*DELL2)*DELI/(2»*DELS) 

SUM = T1*Y(IU) +T?*Y (IUM1) 

IF (SIGMA .EG. 0.) GO TO 9 
CALL SNHCSH (DUMMY*CUl *SIGMAP*DELU1 »2) 

CALL SNHCSH <0UMMY»CU2*SIGMAP*DELU2»2) 

CALL SNHCSH ( DUMMY ♦ CL 1 * SI GM AP*DELL 1 * 2 ) 

CALL SNHCSH (DUMMY*CL2»SIGMAP*DELL2*2) 

CALL SNHCSH ( SS . DUMMY * SI GMAP*DELS »- 1 ) 

SUM = SI|M*TERM(CU1»CL1*T1)*YP(IU) 

* +TERM(CL2»CU2»T2)*YP (IUM1) 

GO TO 10 

9 sum = SUM-Tl * (0FLU2* (DELS+DELUl ) +DELL2* (DELS + DELLl ) ) * 

* YP(IU)/12. 

* -T2* (DELLl* (0ELS*DELL2) ♦DELU1* (DELS+DELU2) ) * 

* YP ( I UMl ) / 1 2 • 


C 



C CORRECT SIGN AND RETURN 
C 

10 CUPVI = SSIGN*SUM 
RETURN 
END 
C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

FUNCTION INTRVL (T,X»N) 

c 

INTEGER N 
REAL T * X ( N ) 

C 

C FROM THE SPLINE UNDER TENSION PACKAGE 

C CODED BY A. K. CLINE AND R. J. RENKA 

C DEPARTMENT OF COMPUTER SCIENCES 

C UNIVERSITY OF TEXAS AT AUSTIN 

C 

C THIS FUNCTION DETERMINES THE INDEX OF THE INTERVAL 
C (DETERMINED BY A GIVEN INCREASING SEQUENCE) IN WHICH 
C A GIVEN VALUE LIES. 

C 

C ON INPUT — 

C 

C T IS THE GIVEN VALUE. 

C 

C X IS A VECTOR OF STRICTLY INCREASING VALUES. 

C 

C AND 
C 

C N IS THE LENGTH OF X (N .GE. 2). 

C 

C ON OUTPUT — 

C 

C INTRVL RETURNS AN INTEGER I SUCH THAT 
C 


c 

I = 1 

IF 



T .LE. 

X (2) ♦ 

c 

I = N-l 

IF 

X (N-l ) 

.LE. 

T 

* 

c 

OTHERWISE 


X (I) 

.LE. 

T .LE. 

X ( I + 1 ) » 


C 

C NONE OF THE INPUT PARAMETERS ARE ALTERED. 
C 

C 

c 

TT = T 

IF (TT .LE. X (2) ) GO TO 4 
IF (TT .GE. X ( M— 1 ) ) GO TO 5 
IL = 2 
TH = N-l 



ooooooooooooooo 


C 

C LINEAR INTERPOLATION 
C 


C 

C 

C 


11= IL+IFIX (FLOAT ( IH-IL) * (TT- 
IF ( TT .LT. X (I) ) GO TO 2 
I p (TT .LE. X ( I + 1 ) ) GO TO 3 


X(IL) )/( 


X ( 


IH)-X(IL) ) 


) 


TOO HIGH 


IL = 1*1 
GO TO 1 
C 

C TOO LOW 
C 

2 IH = I 
GO TO 1 

3 INTRVL = I 
RETURN 

C 

c left end 
c 

4 INTRVL = 1 
RETURN 

C 

C RIGHT END 
C 

5 INTRVL = N-l 
RETURN 

END 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

SUBROUTINE SNHCSH (SINHM, COSHM,X, ISW) 


INTEGER ISW 

REAL SINHM, COSHM,X 


THIS SUBROUTINE 
SINHM(X) 
COSHM(X) 

AND 


UINDER TENSION PACKAGE 
CODED BY A. K, CLINE AND R, J, RENKA 
DEPARTMENT OF COMPUTER SCIENCFS 
UNIVERSITY OF TEXAS AT AUSTIN 


RETURNS APPROXIMATIONS TO 
SINH(X)-X 
COSH ( X ) -1 


COSHMM(X) = COSH (X) -l-x*X /2 
with RELATIVE ERROR LESS THAN 3.42E-1 


4 


ON INPUT — 


o o 


NJ 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


X CONTAINS THF VALUE OF THE INDEPENDENT VARIABLE. 
ISW INDICATES THE FUNCTION DESIRED 


s 

-1 

IF 

ONLY 

SINHM 

IS DESIRED. 

= 

0 

IF 

BOTH 

SINHM 

AND COSHM ARE DESIRED, 

= 

1 

IF 

ONLY 

COSHM 

IS DESIRED, 

s 

2 

IF 

ONLY 

COSHMM 

IS DESIRED, 

ON OUTPUT — 

3 

IF 

BOTH 

SINHM 

AND COSHMM ARE DESIRED. 


SINHM CONTAINS THF VALUE OF SINHM(X) IF ISW .LE. 0 OR 
ISW .EC. 3 (SINHM IS UNALTERED IF ISW .EQ.l OR ISW .EO. 
2 ) . 

COSHM CONTAINS THF VALUE OF COSHM(X) IF ISW .EQ. 0 OR 
ISW .FO. 1 AND CONTAINS THE VALUE OF COSHMM(X) IF ISW 
.GE. 2 (COSHM IS UNALTERED IF ISW .EQ. -1). 


AND 

X AMD ISW ARE UNALTERED. 


DATA SP4/4.50217693381333E-08/, 

* SP3/8. 9527 854421 6390 E-06/. 

* SP2/8, 7204897679 1502E-04/. 

* SP1/4.36314556981690E-02/, 

* SQ1/-6.36854430175110E-03/ 

DATA CP4/1 .78419567490190E-07/, 

* CP3/2.87277229799044E-05/, 

« CP2/2.15151519902028E-03/. 

* CP1/7.58181822756256E-02/, 

* CQ1/-7.51515105679867E-03/ 

DATA 7P3/5.59297116264720E-07/* 

* ZP2/1.77943488030894E-04/, 

* ZP1/1.69800461894792E-02/, 

* ZQ4/1.33412535492375E-09/, 

* 703/-5.80858944138663E-07/, 

* ZQ2/1 .27814964403863E-04/ » 

* ZQ1/-1.63532871439181E-02/ 

XX = X 

AX = ABS(XX) 

XS — X X*X X 

IF ((AX .GE. 2.70) .OR. (AX .GE. 1.15 .AND. 

* ISW .ME. 2)) EXPX = EXP ( AX ) 


SINHM APPROXIMATION 



ooooo ooo ooo 


IF (ISW .EQ. 1 .OR. ISW ,EQ. 2) 60 TO 2 
IF (AX .GE. 1.15) GO TO 1 

SIMHM = ( ( ( ( (SP4*XS+SP3)*XS+SP2)*XS*SP1)*XS+1.)*XS*XX) 

* /{ <SQl*XS+l.)*6.) 

60 TO 2 

1 SIMHM = -( ( (l./FXPX+AX) + AX ) -EXPX ) /2 . 

IF (XX .LT. 0.) SIMHM = -SINHM 

COSHM APPROXIMATION 

2 IF (ISW .NE. 0 .AND. ISW .NE. 1) 60 TO 4 
IF (AX .GE. 1.15) GO TO 3 

COSHM = ( ( ( ( (CP4*XS+CP3) *XS+CP2) *XS+CP1 ) *XS + 1 . ) *XS) 

* /((C01*XS*1.)*2.) 

GO TO 4 

3 COSHM = ( (1./EXPX-2.) +EXPX) /2. 

COSHMM APPROXIMATION 

4 IF (ISW .LE. 1) RETURN 

IF (AX .GE. 2.70) GO TO 5 

COSHM = { ( ( (ZP3*XS*ZP2)*XS+ZPl)*XS+l.)*XS*XS)/( ( ( ( (ZQ4 

* *XS+7Q3)*XS+ZQ2)*XS*ZQ1) *XS*1 . ) *24. ) 

RETURN 

5 COSHM = ( ( (l./EXPX-2.)-XS) ♦EXPXJ/2. 

RETURN 

END 


THE SUBROUTINE AVERAGES THE TOP AND BOTTOM WALL VALUES OF DYJ 

SUBROUTINE A VGDY J ( N J , X J » DY J ) 

DIMENSION XJ (NJ) ,DYJ (NJ) 

N=NJ/2 

IF ( (2*N) .NE.NJ5GO TO 20 
DO 10 1 = 1, N 
SIGN=1.0 

IF (DYJ (I) .LT.0.0) SIGN=-1.0 
DYA=(ARS(DYJ(I) ) +ABS (DYJ ( I+N) ) )/2.0 
DYJ ( I ) =DYA*SIGN 
DYJ ( I +N ) =— DYJ ( I ) 

10 CONTINUE 

W WRITE (6.R9) 

DO 15 1 = 1. N 

WRITE (6.98) I .XJ ( I ) ,DYJ(I) ,XJ(I+N) ,DYJ(I+N) 

15 CONTINUE 
RETURN 

20 WRITE (6.97) 

RFTURN 



o o o 


97 FORMAT ("l"///* RF QUEST FOR AVERAGING DENIED*/ 

$* NJ IS AN ODD NUMBER OF JACKS*/) 

98 FORMAT (2X12,4 (2XG12. 6) ) 

^ 99 FORMAT ("1"//* RESULTS OF AVERAGING DYJ OF TOP AND* 

"R* BOTTOM WALLS*//3X*I*7X*XJ (TOP) *6X*DYJ (TOP) *7X*XJ(BOT) * 

$6X*DYJ (ROT) */) 

END 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

THIS SUBROUTINE WRITES THE PREDICTED RESULTS ON TO TAPE6 

SUBROUTINE NE WY J ( I RL AX ,NRDY , RDY ) 

REAL M,MSET 

INTEGER TEST »RUN, POINT 

COMMON/JACK/N J,NJT »NJR*XJ(100) * Y J ( 1 00 ) » UJ ( 1 0 0 ) , V J ( 1 0 0 ) » 

$U JN (10 0) ♦VJN(IOO) » Y JN ( 1 0 0 ) ,DYJ(100) 

COMMON/P ARAM/HEAD (8) , TEST ♦ RUN , POINT * M * RN * ALPHA » TEMP , RL AX ( 10 ) * 
SNRLAX»PT.PINF»IJINF*IAVG»IS(4) , I0PT(4) »EC(4) , GAMMA ,MSET 
DIMENSION RDY(NPDY) 

DO 10 1=1 »NJ 

RDY ( I ) =RL AX (I RLAX ) *DY J ( I ) 

YJN(I) =YJ(I) +RDY (I) 

10 CONTINUE 

RMSDYTC=RMS (NJT ,RDY ( 1 ) ) 

RMSOYBC=RMS (NJB,RDY (NJT*1 ) ) 

WPITE (6*983) RLAX ( I PL AX) 

WRITE (6*984) (HEAD ( I ) *1=1*8) 

WRITE(6,985) TEST » RUN » POINT 

WRITE (6*986) MSET»RN» ALPHA » TEMP » PT » P I NF , U INF 
WRITE (6*987) 

C WRITE (7* 991) TEST* PUN » POINT »RLAX (I RLAX) »MSET 

00 915 I = 1 * N J 

WRITE (6,988) I * X J ( I ) » Y J ( I ) » Y JN ( I ) * RDY ( I ) 

C WRITE (7,990) XJ (I) ,YJ (I) » RDY ( I ) ,YJN(I) 

915 CONTINUF 

WRITE (6,989) RMSDYTC»RMSDY8C 

983 FORMAT ("1"//* UNSCALED WALL CORRECTIONS FOR A RELAXATION* 

S* FACT0R=*F1 0 ,5) 

984 FORMAT (//3X8A10) * T / \ 

985 FORMAT (3X,*USING RESULTS FROM TEST *,I4,* RUN *,I4,* POINT *14) 

986 FORMAT (//* 3X » ^TUNNEL PARAMETERS ARE** /* 1 0X**M=*j F8#4f 4X**RN= ♦ G9*4 
S,4X,*ALPHA = *,F8.3,3X,*TEMP(DEG P ) =* , F8 . 3 , / , 1 OX , *PT = * ,G9 .4 , 3X , 
S*PINF=**G9,4»1X,*UINF = *»E11.5) 

987 FORMAT (//4X*I*6X*XJ*11X*YJ* 11 X*YJN*1 1X*DYJ*/) 

988 FORMAT (2X13,4 (2XG12.6) ) 

989 FORMAT (//33X*RMS(DYJ(T0P) )=*G12.6/ 

% 33X*RMS (DYJ (BOT) ) =*G12.6) 



990 FORMAT (4F10. 5) 

991 J?nT ( * TEST= * 13 * R UN=*I3* P0INT=*I3* RLAX=*F6.4* MSET=* 
1 0 • 5 ) 

END 


Ul 


APPENDIX C 


SAMPLE CASE 


INPUT DATA LISTING 


6X19 INCH TPANSONIC TUNNEL FLEXIBLE WALL TEST 


39 

3.0728 

31 

35 

0 . 

508.47 

20.616 

12 

12 

19 1 

19 1 

2 0 

2 

• 

1 

2 1 

0 . 

0 . 

0 . 

.25000 


1.0000 




-29.00000 9.50000 

-24.00000 9.49765 

-20.00000 9.49739 

-16.00000 9.50130 

-12.00000 9.50918 

-8.00000 9.51911 

-4.00000 9.55012 

0.00000 9.57337 

4.00000 9.54136 

8.00000 9 .50557 

12.00000 9.49375 

16.00000 9.49179 

-29.00000 -9.50000 

-24.00000 -9.49765 

-20.00000 -9.49739 

-16.00000 -9.50130 

-12.00000 -9.50918 

-8.00000 -9.51911 

-4.00000 -9.55012 

0.00000 -9.57337 

4.00000 -9.54136 

8.00000 -9.50557 

12.00000 -9.49375 

16.00000 -9.49179 

-26.00000 9.50000 -.00900 

-18.00000 9.50000 -.00340 

-14.00000 9.50000 -.00300 

-10.00000 9.50000 -.00830 

-7.00000 9.50000 -.00440 

-6.00000 9.50000 .00110 

-5.00000 9.50000 .00020 

-3.00000 9.50000 -.03590 

-2.00000 9.50000 -.06440 

-1.00000 9.50000 -.09820 

1.00000 9.50000 -.08770 

2.00000 9.50000 -.05520 


914 .76700 

0 0 11 



3.00000 

9.50000 

• • 0 ? 04 0 

5.00000 

9.50000 

•00300 

6.00000 

9.50000 

-.00480 

7.00000 

9.50000 

-•00780 

10.00000 

9.50000 

.00090 

14.00000 

9.50000 

.00610 

18.00000 

9.50000 

.00470 

18.00000 

0.00000 

.00470 

18.00000 

- 9.50000 

.00470 

14.00000 

- 9.50000 

.00610 

10.00000 

- 9.50000 

.00090 

7.00000 

- 9.50000 

-.00780 

6.00000 

- 9.50000 

-.00480 

5.00000 

- 9.50000 

.00300 

3.00000 

- 9.50000 

-.02040 

2.00000 

- 9.50000 

-.05520 

1.00000 

- 9.50000 

-.08770 

- 1.00000 

- 9.50000 

-.09820 

- 2.00000 

- 9.50000 

-.06440 

- 3.00000 

- 9.50000 

-. 03590 , 

- 5.00000 

- 9.50000 

.00020 

- 6.00000 

- 9.50000 

.00110 

- 7.00000 

- 9.50000 

-.00440 

- 10.00000 

- 9.50000 

-.00830 

- 14.00000 

- 9.50000 

-.00300 

- 18.00000 

- 9.50000 

-.00340 

- 26.00000 

- 9.50000 

-.00900 

- 26.00000 

0.00000 

-.00900 

INF = .770943 

UINF 

= 805.542 




NEW JACK PERTURBATION VELOCITY CALCULATIONS (SCALED) 


I 

XJ 

YJ 

UJ 

UJN 

VJ 


VJN 

DVJ 

DU J 

1 

-29,0000 

6.0506 

2.03 

2.03 

-• 2802E- 

14 

-•2802E-14 

0 . 

0. 

2 

-24,0000 

6.0491 

1.18 

-1.99 

-.2428 


- .931 3E-0 1 

.1497 

-3.171 

3 

-20.0000 

6.0489 

.74 

-2.15 

.1943 


. 7926E-0 1 

-.1151 

-2.893 

4 

-16.0000 

6.0514 

.41 

-2.21 

.8707 


.7597 

-.1111 

-2.620 

5 

-12.0000 

6.0564 

.90 

-1.96 

.8623 


1.722 

.8594 

-2.856 

6 

-8.0000 

6.0628 

1.31 

-2.62 

2.542 


2.012 

-.5299 

-3.925 

7 

-4.0000 

6,0825 

2.17 

1.04 

4.754 


8.084 

3.329 

-1.132 

8 

0.0000 

6.0973 

17.21 

5.97 

-.6386 


-1.734 

-1.095 

-11.24 

9 

4.0000 

6.0769 

.02 

.49 

-5.611 


-7.729 

-2.118 

.4667 

10 

8.0000 

6.0541 

1.03 

-3.12 

-3.022 


-3.818 

-.7957 

-4.150 

11 

12.0000 

6.0466 

-.83 

-2.57 

-.6345 


-2.634 

-2.000 

-1.734 

12 

16.0000 

6.0454 

-.98 

-1,93 

.2564 


-1.455 

-1.711 

-.9525 

13 

-29.0000 

-6.0506 

2,03 

2.03 

• 2802E- 

14 

•2802E-14 

0. 

0 . 

14 

-24.0000 

-6.0491 

1.18 

-1.99 

.2428 


•9313E-0 1 

-.1497 

-3.171 

15 

-20 • 0000 

-6.0489 

.74 

-2.15 

-.1943 
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